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PULSED  REFLECTION  AND  TRANSMISSION 
FOR  A  DISPERSIVE  HALF  SPACE 


1.  INTRODUCTION 

This  report  investigates  the  one-dimensional  (1-D)  propagation  of  normally  incident,  pulsed, 
electromagnetic  plane  waves  on  a  half  space  that  is  isotropic,  spatially  homogeneous,  and  dispersive. 
The  dispersive  nature  of  the  medium  is  characterized  by  a  singly  resonant  Lorentz  model  of  the 
refractive  index.  Results  of  previous  studies  of  propagation  at  optical  frequencies  are  extended  to 
include  propagation  at  nonoptical  frequencies  through  media  for  which  a  Lorentz  representation 
is  applicable.  In  the  work  of  Sommerfeld  [1],  Brillouin  [2,3],  and  Oughstun  and  Sherman  [4-8], 
the  signal  is  specified  on  the  medium’s  boundary,  and  propagation  is  analyzed  inside  the  Lorentz 
medium  for  a  constant  transmission  coefficient  of  unity.  Consequently,  these  formulations  account 
neither  for  reflection  from  the  boundary  nor  for  the  frequency  dependence  of  the  transmission 
coefficient.  The  need  for  a  more  general  formulation  that  includes  frequency-dependent  reflection 
and  transmission  coefficients  is  justified  by  rigorous  and  plausibility  arguments.  Not  only  does  the 
approach  delineated  herein  include  reflection,  but  it  also  provides  the  framework  for  addressing 
sources  that  are  external  to  the  medium.  Further,  the  ability  to  determine  the  reflected  field  has 
the  advantage  of  providing  a  diagnostic  capability  for  learning  about  the  medium  through  the 
behavior  of  the  reflected  signal.  The  motivation  for  this  study  stems  from  a  desire  to  apply  the 
uniform  asymptotic  methodology  of  Oughstun  and  Sherman  to  propagation  at  radio  frequencies. 
Specifically,  one  would  like  to  apply  the  uniform  asymptotic  approach  to  two  radar  problems: 
(1)  transionospheric  propagation  at  frequencies  between  10  MHz  and  3  GHz  and  (2)  low-altitude 
propagation  above  the  sea  surface  for  frequencies  between  10  MHz  and  100  GHz.  Of  course, 
the  relations  among  the  physical  parameters  that  describe  these  media  are  different  from  their 
relationships  at  optical  frequencies.  A  byproduct  of  this  analysis  is  its  applicability  to  transmission 
of  short-pulse  signals  by  ultrawideband  (UWB)  radars  [9],  a  recent  area  of  active  research  [10-14]. 

In  Section  2,  the  problem  is  stated,  and  the  reflected  and  transmitted  fields  in  the  frequency 
domain  for  an  arbitrary  pulsed  plane  wave  are  given.  Integral  expressions  of  the  corresponding 
temporal  fields  are  obtained  in  a  straightforward  manner  by  taking  the  inverse  transform  of  the 
frequency-domain  fields;  and  the  geometry,  associated  with  the  transforms  and  with  the  analytical 
model  for  the  refractive  index  n  of  a  singly  resonant  Lorentz  medium,  is  discussed.  In  the  Lorentz 
model  of  the  medium,  the  resonance  frequency  o>o  is  assumed  to  exceed  the  absorption  frequency 
S.  The  remainder  of  the  report  addresses  the  evaluation  of  the  reflected  and  transmitted  fields  for 
an  incident  pulsed  plane  wave  that  is  a  sine  function  with  a  finite  number  of  cycles. 

An  exact  expression  for  the  reflected  field  is  obtained  in  Section  3.  This  expression  is  then 
reformulated  to  provide  more  numerically  tractable  results,  which  are  evaluated  for  Brillouin’s 
choice  of  the  Lorentz  medium  parameters  for  several  values  of  the  sine  function’s  carrier  frequency 
wc.  These  results  are  used  to  obtain  a  physical  understanding  of  the  impact  of  the  dispersive 
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medium  on  the  reflected  field.  In  particular,  one  observes  the  origin  of  the  broadening  of  the 
reflected  field.  Moreover,  the  authors  demonstrate  that  the  reflected  field  is  significant  and  should 
not  be  ignored,  even  at  optical  frequencies  for  certain  ranges  of  uc  and  u>q.  In  addition,  the  early- 
time  and  high-frequency  fields  are  derived  and  compared,  and  the  high-frequency  field  is  shown  to 
be  a  generalization  of  an  earlier  result  by  Colby  [15]. 

In  Section  4,  the  time-domain  transmitted  field  is  treated  under  the  uniform  asymptotic 
methodology  developed  in  Ref.  4  for  evaluating  the  solution  in  the  far  field.  In  this  approach, 
the  observation  points  are  taken  sufficiently  far  into  the  Lorentz  medium  to  guarantee  the  asymp¬ 
totic  representation  of  the  field.  As  mentioned  earlier,  this  analysis  generalizes  optical  formulations 
of  propagation  relative  to  a  Lorentz  medium  to  electromagnetic  propagation  through  that  medium 
by  including  a  frequency-dependent  transmission  coefficient.  Hence  the  extension  is  applicable  in 
situations  where  the  frequency  dependence  of  the  transmission  coefficient  is  important,  such  as 
propagation  at  microwave  and  millimeter-wave  frequencies,  propagation  involving  ultrawideband 
signals,  or  optical  propagation  for  certain  ranges  of  the  physical  parameters.  The  advantage  of 
Oughstun  and  Sherman’s  methodology  is  that  it  partitions  the  inversion  integral  into  three  basic 
components,  which  are  handled  separately.  The  components  roughly  correspond  to  contributions 
from  saddle  points  far  from  the  origin  (Sommerfeld  precursor),  from  saddle  points  close  to  the 
origin  (Brillouin  precursor),  and  from  poles  of  the  inversion  integrand  (main  signal).  The  impact  of 
the  frequency- dependent  transmission  coefficient  is  demonstrated  through  detailed  analysis  of  the 
Sommerfeld  precursor  and  a  heuristic  discussion  of  the  main  signal.  Since  the  objective  here  is  to 
establish  the  credibility  of  including  the  transmission  coefficient  for  propagation  through  a  Lorentz 
medium,  complete  treatments  of  the  Brillouin  precursor  and  the  main  signal  are  not  undertaken. 

Since  the  transmitted  field  depends  crucially  on  the  behavior  of  the  saddle  points  of  the  in¬ 
version  integral’s  generalized  phase,  the  eighth-degree  polynomial  [4,8]  of  which  the  saddle  points 
are  roots  is  corrected.  Although  the  polynomial  has  errors  that  lead  to  incorrect  values  for  the 
saddle  points,  no  errors  are  found  in  any  calculations  involving  the  saddle  points  in  Refs.  4  and  8. 
Moreover,  for  Brillouin’s  values  of  the  physical  parameters,  the  contributions  of  the  distant  saddle 
points  to  the  transmitted  field  are  determined,  the  effect  of  including  the  transmission  coefficient  is 
discussed,  and  aspects  of  the  transmitted  field  are  compared  to  the  reflected  field.  In  particular,  the 
peak  value  of  the  reflected  field  is  shown  to  be  21%  of  the  incident  field’s  amplitude,  and  the  com¬ 
ponent  of  the  transmitted  field  associated  with  the  poles  of  the  inversion  integral  (main  signal)  is 
heuristically  argued  to  be  79%  of  the  field  obtained  when  the  transmission  coefficient  is  unity.  The 
purpose  of  this  example  is  to  demonstrate  that  the  reflection  and  transmission  coefficients  should 
not  be  ignored.  Lastly,  the  Appendix  presents  detailed  investigations  of  the  refractive  index,  re¬ 
flection  coefficient,  and  transmission  coefficient  over  the  complex  plane.  In  addition  to  analytical 
arguments,  two-dimensional  (2-D)  and  three-dimensional  (3-D)  graphs  are  used  to  characterize  the 
behavior  of  each  quantity. 


2.  DERIVATION  OF  FIELDS  IN  TIME  AND  FREQUENCY  DOMAINS 


A  pulsed  plane  wave  E'(z,t),  traveling  in  the  positive  z  direction  in  free  space  (z  <  0),  is 
normally  incident  on  the  dispersive  half  space  (z  >  0)  at  observation  time  t  =  0  (Fig.  1).  The  fields 
in  both  half  spaces  satisfy  Maxwell’s  equations 


dEy  dHx 

-d7  =  ^— 


dHx  dDy 
dz  dt 


dt 


and 


(1) 
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Let  Ey(z,v),  Dy{z, a?),  and  Hx{z,u)  be  the  complex  Fourier  transforms  of  Ey(z,t),  Dy(z,t),  and 
respectively,  for  the  complex  variable  a>.  If  G  represents  any  of  the  temporal  functions 
and  G  is  the  associated  Fourier  transform,  the  complex  Fourier  transform  pair  is 

1  f+°°  f 

<?(u)  =  2tt  dt'  G(*)  =  Jc  (2) 

The  contour  C  is  the  horizontal  line  in  the  complex  a?  plane  that  runs  from  -oo  +  ia  to  +oo  +  ia . 
The  real  constant  a  is  chosen  so  that  all  of  the  integrand’s  singularities  lie  below  the  line  and  so 
that  e  aiG(t)  is  absolutely  integrable  for  all  of  the  fields.  The  existence  of  the  transform  pair  is 
guaranteed  for  any  positive  a  satisfying  these  conditions. 


iilllilll 

— FREE  SPACE  (z<0) 

:  DISPERSIVE  MEDtUM{z>0) 

— ►-  [Fo,  eol 

£(«)  =€o£^(0)  *  £on2(<u)l  . 

I Z 

Fig.  1  -  Geometry  and  specified  coordinate  frame  for  a  pulsed  plane  wave  Ely{z,t)  in  free  space 
(Z  <  0)  that  depends  only  on  its  propagation  direction  z  and  is  normally  incident  on  the 
dispersive  medium  at  Z  =  0.  The  positive  X  direction  points  into  the  paper,  and  the  dispersive 
medium  occupies  the  half  space  Z  >  0.  Free  space  is  characterized  by  the  permittivity  €q  and 
permeability  po,  and  the  dispersive  medium  is  defined  by  the  free-space  permeability  po  and  the 
frequency-dependent  relative  permittivity  €r(u>).  The  quantity  ft  is  the  dispersive  medium's 
refractive  index. 


Transforming  Eqs.  (1)  and  substituting  i)y(z,u>)  =  e(uj)Ey(z,  yield  the  one  dimensional 
Helmholtz  equation  in  the  Fourier  domain  for  the  electric  field, 


d2Ey 
dz 2 


+  fc2  Ey  =  0, 


where  k2{u)  =  fi0€( u)u>2.  The  magnetic  field  is  related  to  the  electric  field  through 


Hr 


i  dEy 
Hqu  dz 


The  permittivity  and  wave  number  are  given  by 


f  «o,  2  <  0, 

1 60€r(w),  z>0, 


(  *o(w), 

\  k(u)  =  k0(u)n(u), 


z  <  0, 
z  >  0, 


(3) 

(4) 

(5) 
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where  c  =  1/^/^ofo  is  the  speed  of  light  in  free  space,  fco(w)  =  w/c  is  the  free-space  wave  number, 
eo  is  the  permittivity  of  free  space,  eT  (=  n 2)  is  the  permittivity  of  the  dispersive  medium,  and  n 
is  the  refractive  index  of  the  dispersive  medium. 


The  analytical  model  for  the  refractive  index  of  a  singly  resonant  Lorentz  medium  is  [16] 

,  s  L  b2  11/2  _  f(w-<4)(« -«/.)] 1/2  (f.s 

L>2-u2  +  2i6u\  ~  [(«-«+ )(«-«_)]  ‘  K  > 

The  quantities  6,  a;o,  and  6  are  positive  real  numbers  and  represent  the  plasma  frequency,  resonance 
frequency,  and  damping  constant  of  the  medium,  respectively.  Since  $  <  u o>  the  branch  points  u>+, 
cj^_,  o;+,  and  lie  in  the  lower  half  of  the  complex  u>-plane  along  the  line,  Im[u;]  =  — <5,  and  are 
given  by 

w+  =  Jul  -  S2  -  iS , 


=  -Vwo  ~£2  ~ 


-  S2  -  iS , 

=  —  yju\  —  £2  —  z5, 

where  Re[u>]  =  u/,  Imfcj]  =  a;",  and  =  u>q  +  i>2  (Fig.  2). 


go" 


Fig.  2  -  Cut  complex  plane  associated  with  the  Lorentz  model  of  the  refractive  index  71  for 
6  <  U o,  where  S  and  U7o  are  the  attenuation  and  resonance  frequencies,  respectively,  of  the  model. 
The  four  vectors,  {u;  —  o/_,u;  —  CJ_,U>  —  U?+,07  —  represent  the  vectors  from  the 

branch  points  {£L7^,U7_,U7_|-,U?!f_}‘  of  71  to  an  arbitrary  point  U)  (=  LJf  +  W,f)  in  the  complex  plane. 
The  magnitudes  and  angles  corresponding  to  these  vectors  are  {// ,rf ,  r,p}  and  {af ,  ,  Tp ,  Ot} > 

respectively.  The  branch  cuts  are  the  horizontal  line  segments  along  U  =  —  iS  that  connect  to 
U7^  and  U7_  to  LO!_. 


In  the  following,  the  superscripts  i,  r,  and  t  are  used  to  denote  the  incident,  reflected,  and 
transmitted  fields,  respectively.  This  investigation  assumes  that  the  incident  field  is  an  arbitrary 
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pulsed  plane  wave  with  leading  edge  that  first  reaches  the  dispersive  medium’s 
at  time  t  =  0.  For  z  <  0  and  t  >  0, 

boundary  {z  —  0) 

Ey(z,u)  =  E'y(z,u )  +  Ery(z,u), 

(8a) 

Hx(z,u)  =  Hi(z,u)  +  Hrx(z,u>), 

(8b) 

Ei(z,u)  =  A(u)eiko^z, 

(8c) 

ETy{z,u)  =  A{u)R{u)e~ik°^z , 

(8d) 

££(*,«)  =  *o(a,)A( u)eik°^\ 

PoL> 

(8e) 

Hrx(z,u)  =  fco(w)  A(u)R(u)e~ik^z; 

HqL} 

(8f) 

and  for  z  >  0  and  t  >  0, 

Ey{z,u)  =  £$(*,«)  =  A(u)T(u)eiko^n^z, 

(9a) 

Hx(z,u)  =  Hl{z,u)  =  -^j^-n(u)A(oj)T(u)eiko^n^z. 

(9b) 

The  reflection  (.R)  and  transmission  (T)  coefficients  are  given  by  the  usual  expressions,  and  the 

function  A  is  the  complex  Fourier  transform  of  £*(0,f).  Thus 

i  r+°°  . 

^(W)  =  2n  J  „ 

(10a) 

R(r  _  1  -  «(«) 

R{U})~  1  +  «(«)• 

(10b) 

r(“)=  i +!(.)• 

(10c) 

Inverting  Eqs.  (8(d)),  (8(f)),  (9(a)),  and  (9(b))  yields 

Ery{z,t)  =  f  A(w)R(u)e-iuWz/c»  du, 

Jc 

(11a) 

Ey(z,t)  =  [  A(u)T(uj)eiuj{z/cKn^-(-ct/zK  fa, 

J  c 

(lib) 

HTx(z,t)  =  —  /  A(u)R(u>)e-iu(t+(z/c»  fa, 

Po  C  Jc 

(11c) 

Hl(z,t)  =  — —  /  n(a;)A(w)T(w)e,'<*’(2/c)ln(w)-(ct/z)]dw. 
PoC  Jc 

(lid) 

These  expressions  differ  from  the  usual  optical  solutions  in  two  respects.  First,  Brillouin,  Sommer- 
feld,  and  Oughstun  specify  the  incident  field  on  the  planar  boundary  and  ignore  the  reflected  field. 
With  the  inclusion  of  the  reflected  fields,  Eqs.  (11)  are  applicable  to  a  wider  range  of  problems, 
that  is,  to  any  electromagnetic  propagation  problems  for  which  the  Lorentz  representation  of  the 
refractive  index  is  valid.  The  second  difference  arises  from  including  the  transmission  coefficient  T 
in  the  pertinent  integrands.  In  optical  formulations,  T  corresponds  to  unity,  its  limiting  value  as 
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u  approaches  infinity.  That  T  can  vary  over  the  range  of  integration  for  nonoptical  signals  may 
substantially  affect  the  value  of  the  integrals  in  Eqs.  (11(b))  and  (11(d)).  The  extent  to  which  the 
inclusion  of  the  transmission  coefficient  affects  the  fields  is  partially  answered.  Moreover,  Eqs.  (11) 
simplify  to  the  optical  solutions  as  the  frequency  of  the  incident  field  approaches  optical  frequencies. 

Because  the  square  root  associated  with  the  index  of  refraction  introduces  four  branch  points 
and  the  exponents  of  the  transmitted  fields  (Eqs.  (11(b))  and  (11(d)))  have  at  least  three  saddle 
points,  the  evaluations  of  the  integrals  in  Eqs.  (11)  are  nontrivial  endeavors.  The  difficulty  in 
determining  the  reflected  fields  is  a  consequence  of  having  to  evaluate  integrals  along  branch  cuts. 
Saddle  points  are  not  issues  for  the  reflected  fields  since  neither  integrand  has  saddle  points.  On 
the  other  hand,  the  integrands  associated  with  the  transmitted  fields  possess  three  or  four  saddle 
points.  The  number  and  locations  of  the  saddle  points  are  not  affected  by  the  inclusion  of  R  and 
T  and  vary  with  position  z  of  the  observation  point  in  the  dispersive  medium  and  time  t  of  the 
observation.  The  spatiotemporal  behavior  of  the  saddle  points  and  their  impact  on  the  asymptotic 
representation  of  the  transmitted  field  are  treated  at  length  in  Refs.  [4, 6, 7, 8]. 


Obtaining  the  time-domain  fields  in  Eqs.  (11)  requires  a  more  detailed  specification  of  the 
incident  field  E{y.  Of  particular  interest  are  modulated  sinusoids,  M(t  -  (z/c))  sin(u;c(t  -  (z/c))),  at 
carrier  frequency  ojc.  The  modulation  function  M  is  zero  for  t  <  z/c.  Rather  than  treat  an  entire 
class  of  modulation  functions,  the  reflected  and  transmitted  fields  are  evaluated  specifically  for  the 
rectangularly  modulated  sine.  Consequently, 


£‘(z,t)  =  sin [wc(*  -  ^)]  U(t  -  -  V(t  -  Z-  -  T0) 


(12) 


The  duration  T0  of  the  modulation  is  a  positive  integer  multiple  of  the  sine  function’s  period 
(T0  =  2n7r/u>c  for  n  =  1,2, . . .),  and  U  is  the  unit  step  function 


for  r  <  0, 
for  t  >  0. 


(13) 


Substituting  Eq.  (12)  into  Eq.  (10a)  and  taking  the  complex  Fourier  transform  yield 


-  ASM. 

±in 
2uc  ’ 


U  ±  ±UC, 
U  =  ±U>C, 


where  . 

S(u>)  = - — 

U  —  U)c  U)  +  U>c 


(14) 

(15) 


3.  EVALUATION  OF  REFLECTED  FIELD 

In  this  section,  the  integral  for  the  reflected  field  ETy  is  investigated,  and  an  exact  representation 
of  E*  is  obtained.  Defining  i  =  t  +  (z/c)  and  substituting  Eq.  (14)  into  Eq.  (11(a))  yield 

Ey(z,t)  =  e(z,t  -  To)  —  e(z,t),  (16) 


'(2’r)  =  s/_ 


oo+ta 
oo  +  ta 


e-,u,TR(u>)S(io)cL;. 


where 


(17) 
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Although  the  function  A  is  analytic  in  the  complex  plane,  Eq.  (14)  permits  Ery  to  be  written  as  the 
difference  of  two  contour  integrals  that  have  simple  poles  at  ±toc  on  the  real  axis.  Observe  that  the 
argument  of  each  integrand’s  exponent  has  no  saddle  points. 


For  r  <  0,  the  contour  Ci  is  closed  from  above  by  the  semicircle  Cfu  (Fig.  3)  to  ensure  that  the 
limit  of  the  integral  along  Cfu  goes  to  zero  as  F  approaches  +oo,  which  is  equivalent  to  invoking 

causality.  Since  the  interior  of  contour  C\  +Cfu  has  no  poles,  the  Cauchy  Integral  Theorem  implies 
that 

e(Z,  t)  =  ^Uxn^  ~  J  e-iuTR(u)S(u)  dco  =  0.  (18) 


CO" 


Fig.  3  Closed  Jordan  contour  C\  +  Cfu,  oriented  in  the  counterclockwise  direction,  for  c(z,  7") 
when  r  <  0  consists  of  the  line  segment  C\  and  the  semicircle  Cfu-  The  line  segment  goes 
from  —  F  -{■  ia,  to  F  ~f~  ia:  and  the  semicircle  has  radius  F  and  center  at  w  =  ia. 


Thus  far,  the  branch  points  of  n(u>)  have  not  played  a  role;  however,  they  are  crucial  in 
determining  e(z,  r)  for  r  >  0,  which  is  now  treated.  The  associated  branch  cuts  are  selected  so  that 
contour  C\  is  closed  from  below,  thereby  guaranteeing  the  convergence  of  the  inverse  transform.  In 
particular,  the  branch  cuts  are  taken  to  be  the  horizontal  line  segments  connecting  u>+  to  to'.  and 
u-  to  u>'_ .  The  branch  is  specified  by  making  the  restrictions 


to  —  to'+  =  peta, 
to  -  to+  = 
to  —  u'_  =  p'etQ 
to  —  to-  =  r'e'^‘ 


—IT  <  a  <  7T, 

0  <  V’  <  27 r, 
-it  <  a'  <  7T 
0  <  <  2tt. 


(19a) 

(19b) 

(19c) 

(19d) 


As  indicated  by  Fig.  2,  the  angles  a,  ip,  a',  and  ip'  are  measured  positively  in  the  counterclockwise 
direction  from  the  positive  w-axis.  Table  1  gives  values  of  these  angles  on  the  line  to  =  -iS.  The 
“+”  and  subscripts  of  locations  A,  B,  C,  D,  and  E  denote  the  top  and  bottom,  respectively,  of 
this  line  over  the  indicated  regions.  Under  these  restrictions,  the  refractive  index  and  the  reflection 
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coefficient  may  be  rewritten  as 


R(“)  = 


EJL 

rr' 


+  i 


l  +  l  +  ffi  +  2ygcos 


(20a) 

(20b) 


where  arg  is  the  argument.  See  the  Appendix  for  detailed  analyses  of  n  and  R.  The  choice  of 
parameters  in  Eqs.  (19)  guarantees  that  n  is  analytic  in  the  cut  plane.  In  particular,  observe  that 
n(A+)  =  n(A_),  n(C+)  =  n(C_),  and  n(£+)  —  n(£L). 


Table  1  -  Values  along  the  line  u  =  —i$  in  the  complex  plane  of  the  argument  of  the  refractive 
index  n  and  of  the  angles  a,  a',  y/,  measured  from  this  line  to  line  segments  connecting  the  point 
u  to  the  branch  points  u;+,  u>+ ,  ,  u;+ ,  respectively.  The  letters  A,  £,  C,  D,  E  indicate  the  disjoint 

line  segments  that  comprise  this  line  as  illustrated  in  Fig.  2,  and  the  +  (-)  subscripts  denote  the 
top  (bottom)  of  each  segment. 


Angle  Value 

EH 

A _ 

B+ 

B _ 

C- 

D _ 

£+ 

a 

0 

0 

7 r 

—  7T 

7 r 

—  7T 

7 T 

—  7T 

7T 

—  7T 

i> 

0 

27T 

0 

2;r 

7T 

7T 

7T 

7T 

7T 

7T 

a' 

0 

0 

0 

0 

0 

0 

0 

0 

7T 

—  7T 

0 

27T 

0 

27 r 

0 

27T 

7T 

7T 

7T 

7T 

argfra] 

0 

-27T 

7r/2 

— 5tt/2 

0 

-2tt 

— 7T /2 

-3x/2 

0 

—  27T 

The  contour  C\  is  converted  to  a  Jordan  contour  Cj  by  adjoining  several  contours  (Fig.  4): 
Ct  =  Cfl  +  Y?n= i  Cn-  Applying  Cauchy’s  residue  theorem  and  letting  F  — >  +oo  and  e  0+ 
simultaneously  give 

e(z,r)  =  -27ri[Res(-u>c)  +  Res(u;c)]  -  —  [ Ib+  +  Ib _  +  Id+  +  Id_],  (21) 

where  Res(u>p)  is  the  residue  of  the  integrand  at  uv  and 


Jb+  =  lirn 

+  c-+04-  , 

1  R(u)e~iuTS(u)du, 

(22a) 

Ib  —  lim 

[  R{u>)e~iu,TS(u)du, 

'C3 

(22b) 

Ip,  =  lim 

+  c— ►O-f  , 

f  R{u;)e-iuTS(u)du, 

Jc9 

(22c) 

Id  =  lim 

c-»0+  t 

f  R(uj)e~,u,TS(u)duj. 

/c„ 

(22d) 
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CD" 


Fig.  4  Closed  Jordan  contour  Cfl  +  Si'Ll  Cn,  oriented  in  the  clockwise  direction,  for  e(z,  r) 
when  T  >  0  consists  of  five  horizontal  line  segments  {CUC3,  C5,C9,Cn},  four  vertical  line 
segments  {C2,C7,CS,Ci3},  four  circles  {C^C^CicCia}  of  radius  £,  and  one  semicircle 
Cfl  of  radius  F .  The  four  circles  are  centered  about  the  branch  points  {w'_,  W_,o;+, 

the  horizontal  line  segments  are  the  tops  and  bottoms  of  the  two  branch  cuts,  and  the  Semicircle 
is  centered  at  a?  =  ia. 


To  obtain  an  expression  for  e(z,  r),  the  four  line  integrals  in  Eqs.  (22)  must  be  evaluated,  and 
the  two  residues  must  be  obtained.  Before  evaluating  the  line  integrals,  define  two  more  quantities: 

u>2  =  ful  -  62  and  w3  =  -  62.  (23) 

Since  u\  =  +  b2,  w3  >  w2. 


Along  the  top  of  the  right  cut,  which  corresponds  to  the  region  B+  in  Fig.  2,  parameterize  C5 

as:  cj  -  r  +  w+  =  r  +  w2  -  iS  with  magnitude  r  running  from  u>3  -  w2  -  £  to  £.  The  remaining 
magnitudes  of  Eqs.  (19)  are 


p  =  -  u>2  -  r, 

r'  =  r  +  2u>2, 

pf  =  u?3  +  ^2  +  V. 

From  Table  1,  a  =  tt  and  =  a'  =  V  =  0.  Consequently,  Eqs.  (20)  become 


(24) 


n(u)  =  i  fr*  +U}2  +  r)(a;3  -bfj-  r)l 1/2 

L  r(r  +  2oj2)  J  ’ 

_  V/r(r  +  ^2)  ~  +  k>2  +  r)(o;3  —  q>2  —  r) 

\fr{r  +  2u>2  )  +  i  yj  (w3  +  u;2  +  r)(w3  —  w2  —  r) 


(25a) 

(25b) 


After  substituting  Eq.  (25(b))  into  Eq.  (22(a))  and  taking  the  limit,  one  obtains  Eq.  (26(a))  as  the 
expression  for  IB+.  The  integral  for  IB_  is  gotten  by  the  same  parameterization.  On  the  other 
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hand,  along  the  left  cut  use  the  parameterization:  =  —  r'  +  u;_  =  —  r*  —  u>2  —  iS  with  r9  running 
from  £  to  CJ3  —  Ct?2  £•  Consequently, 


IB+  =  -e-^+i^T  [ 

Jo 

IB_  =  e-<*+iwa>T  / 
Jo 

ID+  =  e-^-iu,2)T  / 

Jo 

ID_  =  - e 


W3— U>2 


1 


1 


U>3  —  U)2 


W3  — 0)2 


[r  +  u>2  —  r  +  o^2  +  o;c  — 

1  1 


R(u )  dr, 


[r  +  u;2  —  wc  -  r  +  u;2  +  uc  -  i6 
1  1 


|_r  +  u;2  +  uc  +  i6  r  +  u>2  -  uc  +  J 


1 14/3  —  U>2 

/  P~irT 

1 

1 

f  e 

'o 

r  +  u?2  +  uc  +  ib 

r  +  u>2  —  vc  +  ib. 

R*(u>)  dr, 
R'(u)dr, 
R(u)  dr. 


(26a) 

(26b) 

(26c) 

(26d) 


Observe  that  /B+  =  [/b+]*  and  Id_  =  [Ib-]*-  Therefore,  IB+  +  IB_  +  Ip+  +  Id-  =  2Re{/B+  + 
IB_  },  so  that  the  field  is  expressed  in  terms  of  integrals  along  the  right  branch  cut  only.  After 
some  algebra, 


{o-  yu/3— U>2 

l.e-V+i-»)T  J 


r  +  u>2  + 


[r  +  w2  —  uc  —  iS 
~  \J r(r  +  2o;2)  (b2  -  r(r  +  2w2 ))  dr. | 


(27) 


It  remains  to  calculate  the  residues  at  ±uc.  First, 

Res(u>c)  =  ±-R(  uc)e-iu‘\  (28) 

47T 

Since  n(— w)  =  [n(u>*)]*,  R(-uc)  =  [R(w*)]*,  which  implies  that  Res(-o>c)  =  ~[Res(wc)]*.  Com¬ 
bining  this  result  with  Eq.  (28)  yields 

— 27ri[Res(— u>c)  +  Res(u;c)]  =  cos(wcr)Im[R(a;c)]  —  sin(u;cr)Re[R(a;c)]  (29) 

as  the  contribution  of  the  poles  to  e(z,r). 


Equation  (18)  and  the  substitution  of  Eqs.  (27)  and  (29)  into  Eq.  (21)  imply  that 
e(z,  r)  =  f  cos(u>cr)Im[R(u>c)]  —  sin(o;£.r)Re[R(a>c)] 


-  ^Hie~,,+i“,)T  L 


UJ3  —  U/2 


r  +  cj>2  —  cac  —  ib  r  +  u  2  +  —  iS  J 


(30) 


X  yjr{r  +  2 w2)(&2  -  r(r  +  2u>2))  dr|^  t/(r). 


To  perform  numerical  computations,  it  is  desirable  to  recast  this  expression.  After  more  algebra, 
one  obtains 

e(z,r )  =  ^cos(u;cr)Im[R(a;c)]  -  sin(a>cr)Re[.R(u>c)] 


—  Sr 


■zb2 


1  \ 

sin(o;2r)/1(r)  -  cos(w2r)/2(r)  U(r), 

J  / 


(31) 
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3-W2  f  1/2 

7i(r)  =  J  ^i(r,r)[r(r  +  2a;2)(62 -r(r  +  2w2))j  dr,  (32a) 

/a'3-a'2  r  ,  i/2 

Mr)  =  yo  Nz(r,T)  |r(r  -f-  2w2)(62  —  r(r  +  2w2))J  dr,  (32b) 

jy  (r  r)  _  2a)4(r  +  ~  fl  ~  ^2]  cos(rr)  +  4uJ(r  +  a;2)  sin(rr) 

[(r  +  a>2  —  uc )2  +  £2]  [(r  +  a;2  +  «c)2  +  £2]  ’  (32c) 

jV2(r  r)  _  ^(r  +  a>2)  cos(rr)  -  2a;e[(r  +  aj2)2  -  a>2  -  6*]  sin(rr) 

[(r  +  -  ucf  +  £2]  [(r  +  u>2  +  wc)2  +  62\  '  (  2  ) 

Evaluation  of  the  real  and  imaginary  parts  of  R(uc)  are  obtained  from  substituting  the  following 
expressions  into  Eq.  (20(b)): 

Afp?\  _  \{ul  -^i)2  +4w2^2l1/4  ,  , 


(33b) 


'Lc*  LK2-^o2)2  +  4^2J  ’ 

f  rr- aretan  0  <  «c  <  w3 


1  2  ’ 


U)c  =  C^3 
^3  < 


(  aretan  >  <*>3  <  wc, 

a'l  =  aretan  f - - — V 

,w_w#  Vw3+a;cy’ 

'  7T- aretan  (^-),  0  <  wc  <  w2 

^L=wc  =  *  J’  OJc  =  CJ2 

.  arctan(;rf^)>  w2  <  u>c, 


aretan 


*'LW.  =  Cretan  ( 


Finally,  substituting  Eq.  (31)  in  Eq.  (16)  yields  the  following  computationally  useful,  exact  expres¬ 
sion  for  the  reflected  field, 

ETy(z,t)  =  \R(uc)\sm(uci -  arg[£(u>c)])  -  U(t-T0)] 

+  {e~5t/(*)  sinM  “  "(OMO 

-  *-«»-*»  j(f  _  To)  sin(w2  (i -  T0)  -  u{t  -  T0))U(i  -  T0)  J,  (34) 


where 


7(r)  —  \/[7i(r)]2  +  [^(r-)]2  and  v(t)  =  aretan  ( -ffl  ^ 


Note  that  the  field  in  the  right  side  of  Eq.  (34)  consists  of  contributions  from  the  poles  (first 
term)  and  the  branch  cuts  (second  term).  The  contribution  from  the  poles  has  duration  equal 
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to  the  incident  pulse’s  duration,  whereas  the  contribution  from  the  branch  cuts  lasts  indefinitely, 
albeit  they  exponentially  dampen  out.  At  a  point  z  <  0,  one  does  not  observe  the  signal  until 
the  observation  time  t  exceeds  —z/c,  which  is  the  time  it  takes  the  signal  to  reach  the  observation 
point  after  reflection  from  the  boundary  (z  =  0)  of  the  dispersive  medium.  Then  it  is  observed 
for  a  time  To,  after  which  it  rapidly  diminishes  to  a  negligible  value  for  t  >  To  +  (—z/c).  As  one 
may  deduce  from  the  second  term  in  Eq.  (34)  and  the  derivation  preceding  it,  although  the  field 
diminishes  exponentially  as  and  it  does  not  vanish  identically  for  t  >  To  +  (-z/c). 

The  contribution  for  t  >  To  +  (—z/c)  arises  only  from  the  integrals  along  the  branch  cuts.  If 
no  branches  were  present,  the  reflected  field  would  have  the  same  duration  as  the  incident  field. 
Consequently,  it  is  the  presence  of  the  branch-cut  integrals  that  is  directly  responsible  for  broadening 
the  time  extent  of  the  reflected  field.  The  decay  rate  of  the  branch-cut  contribution  increases  as 
the  absorption  coefficient  5  increases,  and  the  magnitude  decreases  as  the  square  of  the  plasma 
frequency  b  increases. 

The  periodicity  (To  =  2n7r/u;c)  of  the  sine  function  allows  one  to  express  the  pole  contribution 
in  terms  of  the  window  function  over  [0,To].  Further  note  that  the  pole  contributions  result  in 
a  replica  of  the  incident  signal,  which  has  a  delay  in  phase  as  a  result  of  the  reflection  from 
the  medium’s  boundary  and  which  has  a  magnitude  that  depends  on  the  behavior  of  R  as  a 
function  of  uc  relative  to  6.  According  to  Eq.  (20b),  the  phase  delay,  the  argument  of  the  reflection 
coefficient,  is  negligible  when  the  imaginary  part  of  R(uc)  is  nearly  zero  or  equivalently  when 
sin[(a  +  a'  —  $  -  V>')/2]  =  0.  On  consulting  the  figures  in  the  next  subsection,  one  can  deduce 
that  this  occurs  in  at  least  two  situations,  when  uc  <  min{u;o,£, b}  or  when  uc  >  max{u>o,£, b}.  If 
u?c  has  the  same  order  of  magnitude  as  the  dispersive  medium’s  parameters,  Im[i2(u;c)]  cannot  be 
approximated  by  zero. 

Similarly,  the  two  components  of  the  branch-cut  contributions  are  expressed  in  Eq.  (34)  as 
exponentially  damped  sinusoids  with  nonzero  time-dependent  phases.  Although  the  contributions 
from  the  branch  cuts  are  expressible  individually  as  modulated  sines,  the  two  step  functions  cannot 
be  combined  to  get  a  window  function  (as  the  pole  contributions  are)  because  2tt/u;2  is  not  neces¬ 
sarily  an  integer  multiple  of  To-  The  first  exponential  term  depends  on  i  and  most  affects  the  field 
just  after  t  exceeds  -z/c,  whereas  the  impact  of  the  second  exponential  term  is  not  felt  until  t  is 
greater  than  or  equal  to  To  +  (—z/c).  The  preceding  conclusions  are  illustrated  by  the  following 
example. 


3.1  Example  of  Reflected  Field 

In  this  section,  the  time  histories  of  the  reflected  field,  corresponding  to  eight  choices  of  the 
incident  field’s  carrier  frequency  uc  at  observation  point  z  =  — 10~5  m  are  plotted  for  Brillouin’s 
choice  [17]  of  the  Lorentz  medium  parameters:  u?o  =  4.0  X  1016s_1,  b2  =  20.0  X  1032s~2,  and 
S  =  0.28  X  1016s-1.  Moreover,  the  incident  field  is  taken  as  a  single  cycle  of  the  sine  function;  that 
is,  E'y  has  duration  2tt/ujc  ( n  =  1).  Note  that  this  signal  is  a  short-pulse  ultrawideband  signal  [18]. 
As  Eq.  (34)  and  the  discussion  following  it  indicate,  Ey  depends  on  the  value  of  uc  relative  to  the 
parametric  frequencies  u;o,  6,  and  5.  In  particular,  the  behavior  of  i?(u>c),  as  uc  varies  along  the  real 
line,  is  crucial  to  the  behavior  of  Ey.  Consequently,  the  behavior  of  R  versus  uc  is  briefly  discussed 
before  addressing  Ey. 

Figure  5  plots  various  aspects  of  R(uc).  The  horizontal  axis  is  uc  divided  by  cjq;  thus  unity 
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corresponds  to  uc  =  u?0.  To  gain  additional  perspective,  note  that  =  6  X  1016s“1,  u?2  = 
3.9902  X  1016s-1,  w3  =  5.9935  X  1016s_1,  u±  =  (±3.9902  -  0.28i)  X  10V1,  and  =  (±5.9935  - 
0.28i)  x  1016s_1,  which  are  associated  respectively  with  the  complex  numbers  1.5,  0.9975,  1.4984, 
±0.9975  —  0.07i,  and  ±1.4984  -  0.07i  in  the  figures.  The  phase  of  R(uc )  is  —180°  in  the  limit  as 
UC  -*  0±,  is  strictly  monotonically  increasing,  and  asymptotically  approaches  0°  as  wc  increases 
above  2w0.  In  fact,  for  ojc  in  (0,0.8wo]  or  [1.8o>o,oo),  |arg[^(wc)]|  <  10°.  Consequently,  the  small- 
angle  approximations  to  the  sine  and  cosine  of  the  argument  of  R{uic)  are  fairly  accurate  for  this 
range  of  u>c.  The  magnitude  of  R(uc)  is  0.2  at  uc  =  0,  rises  to  its  peak  at  uc  =  1.3w0,  and  falls 
rapidly  to  0.  Therefore,  the  magnitude  of  the  contributions  from  the  poles  of  the  reflected  field  will 
be  greatest  for  uc  near  1.3u?o.  Although  the  pole  contributions  are  maximized  at  u  =  1.3wo,  the 
reflected  signal  with  greatest  energy  appears  at  the  resonance  of  the  medium. 


Fig.  5  -  Behavior  of  the  reflection  coefficient  R  associated  with  a  single-resonance  Lorentz  medium  vs 
carrier  frequency  U)c:  (a)  real  (light  curve)  and  imaginary  (dark  curve)  parts  of  R,  (b)  magnitude  of  R, 
and  (c)  argument  of  R.  The  Lorentz  medium  is  defined  by  u>0  =  4.0  X  1016s_1,  b2  =  20.0  X  1032s~2, 
and  8  =  0.28  X  1016s-1. 


To  observe  the  distortions  to  the  incident  field  that  are  caused  by  the  dispersive  half  space,  Er 
is  compared  to  the  reflected  field  that  is  obtained  when  the  boundary  is  perfectly  conducting!  In 
that  case,  R(u)  is  identically  equal  to  -1,  the  transmitted  field  is  zero,  and  the  expression  for  the 


14 


Mokole  and  Samaddar 


reflected  field  is 

~[U{i)  -  U(t-T0)]  sin (uj),  (36) 

which  is  what  one  intuitively  expects.  Equation  (36)  and  Ey  are  compared  graphically  in  Fig.  6 
for  eight  values  of  u>c,  where  the  horizontal  axes  are  normalized  time  6  =  ct/\z\,  a  dimensionless 
quantity.  Since  the  reflected  fields  are  not  observed  until  t  >  |z|/c,  nonzero  values  do  not  occur  until 
0  >  1.  The  choices  of  u;c,  the  associated  evaluations  of  the  reflection  coefficient,  and  other  pertinent 
parameters  are  provided  in  Table  2.  For  example,  the  first  column  of  Table  2  permits  one  to  identify 
uc  with  the  abscissa  of  Fig.  6  and  the  value  of  R(wc)-  Also,  since  the  reflected  field  for  the  perfectly 
conducting  boundary  is  nonzero  when  1  <  0  <  1  +  cTol^l"1,  columns  three  and  four  give  the 
duration  and  right  endpoint,  respectively,  of  this  field.  The  smallest  values  (0.0001u?o,0.01a;o,0.1a;o) 
are  chosen  to  represent  the  behavior  of  Ey  for  u ;c<w0,  while  0.25ic?o  *s  chosen  because  it  is  a  value 
that  Oughstun  and  Sherman  use  in  much  of  their  work.  The  next  largest  value  (u>o)  is  chosen  to 
see  the  impact  at  the  dispersive  medium’s  resonance  frequency.  The  sixth  (1.26366cjo)  and  seventh 
(yi3/8a;o)  values  of  uc  are  of  interest  because  they  are  the  locations  of  the  minimum  of  Im [R(uc)] 
and  the  zero  of  Re[i?(u;c)],  respectively.  The  remaining  choice  of  uc  is  selected  to  represent  the 
behavior  of  Ey  for  large  values  of  the  carrier  frequency. 


Table  2  -  Choices  of  angular  carrier  frequency  uc  of  the  single-cycle-sine  incident  field 
for  which  the  reflected  field  Ey  is  plotted  in  Fig.  7,  as  well  as  related  evaluations. 
Column  1  identifies  the  normalized  frequency  uc/uo  with  wc;  columns  3  and  4,  resp., 
are  the  normalized  duration  and  the  normalized  right  endpoint  of  the  support  of  Ey 
when  the  boundary  is  perfectly  conducting;  and  columns  5  and  6  are  evaluations  of 
the  real  and  imaginary  parts  of  the  reflection  coefficient  R  at  u;c. 


Wc/Wo 

<*>C 

(rad/s) 

m 

1  + (cTo/M) 

Re[i?(wc)] 

Im  [R(uc)] 

0.0001 

4  X  1012 

47.09129 

48.09129 

-0.200000 

-0.000002 

0  .01 

4  x  1014 

0.47091 

1.47091 

-0.200013 

-0.000187 

0  .1 

4  x  1015 

0.04709 

1.04709 

-0.201324 

-0.001893 

0  .25 

1016 

0.01884 

1.01884 

-0.208582 

-0.005098 

1 

4  x  1016 

0.00471 

1.00471 

-0.552437 

-0.276488 

1.26366 

5.05464  x  1016 

0.00373 

1.00373 

-0.032778 

-0.755325 

(13/8)1/2 

5.09902  x  1016 

0.00369 

1.00369 

0.000000 

-0.754425 

10 

4  x  1017 

0.00047 

1.00047 

0.003176 

-0.000045 

Normalized  Time 


Normalized  Time 


(C) 


(d) 


(e) 


ff> 


Fig.  6  -  Comparison  of  the  reflected  field  when  the  medium’s  boundary  is  perfectly  conducting  (light  curve)  and  the 
reflected  field  from  the  Lorentz  medium  at  z  =  — 10~5  m  with  u o  =  4.0  X  1016s“1,  62  =  20.0  X  1032s”2, 
and  6  =  0.28  X  1016s  1  (dark  curve)  for  eight  values  of  the  carrier  frequency  Wc:  (a)  O.OOOlc^o,  (b)  O.OluJo,  (c) 
0.1ct>o>  (d)  0.25u>o,  (e)  (f)  1.26366u?o,  (g)  ^/l3/8u;o,  and  (h)  10u>o*  The  horizontal  axes  are  normalized  time 

9  =  ctf\z\. 
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After  studying  Fig.  6,  Table  2,  and  Eqs.  (34)  and  (36),  one  deduces  several  facts  about  the 
overall  behavior  of  Ey  as  uc  varies.  Initially,  for  very  small  wc,  Ey  is  nearly  a  replica  of  the 
undistorted  reflected  field,  scaled  by  a  factor  of  0.2.  Generally,  as  uc  increases,  the  duration  To  of 
the  undistorted  signal  decreases,  and  the  distorted  field  steadily  broadens,  so  that  eventually  it  is 
many  times  wider  than  the  undistorted  pulse.  Concurrently,  the  maximum  value  of  the  distorted 
pulse  increases  from  0.2  to  a  peak  of  roughly  0.5  at  the  medium’s  resonance  frequency  (u>c  =  u;0)  and 
then  decreases  to  zero.  As  the  carrier  frequency  increases,  one  also  observes  leading-edge  distortion 
and  damped  oscillatory  behavior  in  the  field’s  trailing  edge.  Further,  for  very  large  wc,  a  steep 
leading  edge  is  followed  by  many  zero  crossings. 

More  specifically,  for  very  small  values  of  u?c  (0.0001a;o,0.01a;o,0.1u;o),  the  pulse  widths  of  Ey 
and  the  undistorted  field  are  nearly  the  same,  except  for  some  negligible  ringing  in  the  trailing 
edge.  As  uc  increases  to  O.Icjq,  pulse  broadening  and  trailing-edge  ripple  are  more  pronounced, 
and  leading-edge  distortion  becomes  discernible.  All  three  effects  increase  as  uc  approaches  the 
medium’s  resonance  frequency  (cjq),  where  the  signal  attains  its  largest  excursion  (0.5)  from  zero. 
Moreover,  Ey  consists  roughly  of  three  cycles  with  a  duration  of  approximately  thrice  the  undis¬ 
torted  signal’s  pulse  width.  As  uc  increases  above  uo,  the  pulse  broadening  and  trailing-edge  ringing 
are  greater,  but  the  maximum  value  decreases  to  zero.  Even  though  1.26366u>o  and  ^/l3/8uo  cor¬ 
respond  to  special  values  of  relative  to  i?(u;c),  the  behaviors  of  Ey  for  these  values  are  similar  to 
that  for  u>0.  The  behavior  of  Ey  for  large  u>c  is  typified  by  uc  =  IOcjq*  Since  the  un distorted  pulse 
for  this  value  of  the  carrier  has  an  amplitude  that  exceeds  the  peak  of  Ey  by  at  least  fifty  times 
and  occurs  within  the  first  “half  cycle”  of  Ey  (Table  2),  the  undistorted  signal  is  not  plotted. 

Lastly,  for  u?c  <  3u;0,  the  maximum  value  of  Ey  is  a  significant  fraction  (exceeds  0.2)  of  the 
amplitude  of  the  undistorted  reflected  field.  Clearly,  a  nontrivial  amount  of  the  incident  field’s  en¬ 
ergy  resides  in  the  reflected  field  for  these  values  of  uc.  Hence,  the  reflected  field  should  be  included 
in  any  Lorentz  formulation  where  the  carrier  frequency  is  less  than  three  times  the  medium’s  res¬ 
onance  frequency,  regardless  of  what  portion  of  the  electromagnetic  spectrum  is  being  considered. 
In  particular,  Oughstun’s  choice  (0.25gjo)  of  uic  falls  in  this  range  of  frequencies.  Moreover,  because 
the  reflected  field  is  not  negligible,  it  may  be  used  as  a  diagnostic  for  the  dispersive  medium. 


3.2  Early-Time  Field  Behavior 


The  early- time  behavior  of  the  reflected  field  at  observation  point  —  z  corresponds  to  t  slightly 
greater  than  —z/c.  Consequently,  the  reflected  field  is  analyzed  for  —  z/c  <  t  <  —z/c  +  te  (0  < 
t  <  te ),  where  te  is  a  small  positive  number  less  than  T0.  Hence  in  Eq.  (34),  U(t  -  To)  =  0  for 
0  <  i  <  te,  and  the  early-time  reflected  field  is 


££,et(-M)  =  sin(wcf)Re[i2(wc)]  -  cos(wci)lm[R(wc)]j 

sin(w2<)Xi  (t)  -  cos  (cj2t)l2  (t)  j . 


+  ^e 


(37) 


Alternatively,  when  closed-form  solutions  are  unavailable,  one  resorts  to  high-frequency  approx¬ 
imations  of  various  quantities  to  obtain  analytical  expressions.  The  remainder  of  the  subsection  is 
spent  on  deriving  a  high-frequency  representation  of  Ey.  Begin  with  Eq.  (16)  and  recast  it  as 

/oo+ta  _  ^  r  co-Ha 

Q(u)R{u)A(u  -  u>c)e~tuli  du  -  /  Q(u)R(u)A(u  +  u>c)e~,ui  du>,  (38) 

-oo+ia  J  —  oo-f*a 
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where 

@(“0  =  -  *)  Md  A(o>)  =  i.  (39) 

Because  ]n(w*)]*  =  n(-w),  it  follows  from  Eqs.  (10b)  and  (39)  that  R*(u*)  =  R(-u),  A(-u*  - 

Uc\  =  ~A*{u  -  u;c),  and  Q*{u*)  =  Q(-u>).  In  the  second  integral  of  Eq.  (38),  make  the  change  of 
variable  rj  =  o>*  and  apply  the  preceding  properties  to  get 


*?(*.«)  =  Ref  r  “  R(") 

.*/— oo+ia  27T(u;  0>c) 


_  /'00+,a  .R(a; 

J—oo+ia  27r(w  - 


) 


-iu/t 


oo-fia  27t(cj  a>c) 


da? 


(40) 


The  reflected  field  is  now  ready  for  approximate  analysis.  The  Laurent  series, 

-  [J  -  j  0 + B + B) + 

■R(cj)  b2  [q2  a2wc  a4  +  a2u>2  wc(a4  +  a2w2 )  „  /  _7s 

a>-u,c  "  4  [Z?  +  ~  +  — ^ - +  Rfi(u  7) 


.b2 
1  4 


j_!*Vc  ,  as  +  a3 u2  ,  wc(a5  +  a3w2)  ,  „  , 

w4  a;5  +  w6  +  ji  +  "«(w  ) 


are  obtained  after  several  applications  of  the  binomial  series  for  large  u>,  where 

<Pi  =  1, 

<P2  =  u>l  -  AS2  +  ^b2, 

<P3  =  Uq  —  12 S2u%  +  h2u>S  -  3b2S2  +  h4  +  16<54, 

z  o 

ti  =  6, 

i>2 


=  26  (u2  -  262  +  \b2^j  , 

1^3  =  g  —  32 62Vq  +  362o;q  +  — 64^  ’ 


a2  =  v?i, 

a3  =  2^i, 


62 


<*4  =  V>2  +  —  V>i, 

a5  =  2V>2  +  629»i^i, 

=  ^3  +  (  2<£i  <p2  +  —  <fl  -  4^2  J  , 


2<Pl  <P2  + 
a7  =  2-^3  +  &2  (^P\i>2  + 


(41) 


(42) 


(43) 


The  functions  Rn\,  Rno ,  Rfi,  and  Rf2  represent  the  remainders  of  series,  and  their  arguments 
indicate  the  lowest-order  terms  in  the  remainders.  The  high-frequency  representation  for  n(w),  as 
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a  Laurent  series,  effectively  eliminates  the  branch  points  as  singularities.  As  will  be  seen  shortly, 
this  greatly  simplifies  evaluation  of  the  integrals  in  Eq.  (40). 

Since  each  integral  in  Eq.  (40)  has  the  same  general  form,  consider 


/ 


oo+ta 


R(u) 


—oo+ia  2*-(w  -  Uc) 


du. 


Recall  from  the  discussion  regarding  Eq.  (18)  that  C\  is  the  finite  version  of  the  path  C  in  the 
preceding  integral.  By  closing  C\  with  a  semicircle  of  radius  F  from  above  and  by  letting  F  — >  oo, 
the  path  integral  is  zero  for  r  <  0.  For  positive  r,  the  Jordan  contour  is  simply  C\  +  Cfl,  where 
Cfl  is  the  semicircle  of  radius  F  that  connects  the  ends  of  C\  from  below  (Fig.  7).  According  to 
the  Cauchy  Integral  Formula,19 


lim  / 

F->+co  JCl+CFL 


2  7T 


du  =  ^-(-i)n+1rn 
nl 


Substituting  the  high-frequency  representation  of  R(u)(u  -  uc)~l  given  by  Eq.  (42),  integrating 
term  by  term,  and  applying  the  Cauchy  Integral  Formulas  yield 


/ 


oo-f 


R{u) 


b2 


,  — - ,e~iuT  du  =  U{r)- 

-oo+ia  2ff(w  -We)  4 


a2uc  3  a3uc  4  Uc(a4  +  a2ul)  5  uc(a5  +  a3u2)  6 

3!  4!  T  5!  T  +  6!  T 

.2 


+  Ri(r7) 


+  i 


a2.2  “3,3  04  +  a2U2c  _4  a5  +  a3U2  _5  /_6\ 

2!  "3!r - 4!  T  + - 5! - T  +  R^T  ) 


(44) 


where  Ri  and  R2  are  the  series’  remainders  for  the  real  and  imaginary  parts,  respectively. 


Fig.  7  -  Closed  Jordan  contour  C\  +  Cfl,  oriented  in  the  clockwise 
direction,  consists  of  the  line  segment  C 1  and  the  semicircle  CfL' 
The  line  segment  goes  from  —  F  +  ia  to  Jp  +  ia,  and  the  semicircle 
has  radius  F  and  center  at  U)  =  ia. 
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Setting  r  equal  to  i  To  and  i  in  Eq.  (44),  substituting  those  expressions  into  Eq.  (40),  and 
using  Eqs.  (43)  lead  to 


‘>  =  m  -  To)j  [§(<  -  Tof  -  *g(f  -To)*-  |  (ul  -U*  +  ^+  oi)  (i  -  nf 


■  u(i) 


b2 


■  ^  (2a ;2  -  4S2  +  b2+  a,2)  (i  -Tof  +  Ri  (ft  -  Tof) 

COc&  2 SwC24  (  2  AC<1  6 2  9\  ^ 

3!*  + 


+  (2w l  -  4  S2  +  b2  +u2c)F  +  R1  (f  ) 


(45) 


a  high-frequency  representation  of  ETy.  This  result  is  obtained  by  assuming  that  \u\  is  large  enough 
to  ensure  the  convergence  of  the  series  used  to  obtain  R{u){u  —  u>c)— 2.  At  the  very  least,  |w|  > 
max{o;c,a>o,<5,6}.  Moreover,  when  |w|  is  large  enough,  0  <  t  <  te  and  U(t  -  T0)  =  0  Thus  for 
0  <t<te,  ' 


jpr  (  +\  _  ^  2 Su>c?4  uc  f  o  a  c2  ^  9\ 

Ey,hfz^)  ~  4  [  a; 1  ~  — *  ~  (^o  ~  M  +  “cj  t 

+  (2wo  ~  4£2  +  62  +  o>2)  f6  +  .Ri  (f)j . 
In  the  limit  as  the  attenuation  factor  £  goes  to  zero,  Eq.  (46)  reduces  to 

JWM)  =  7  [~lp  + 1  ("» +  t  +  <4)  ?  +  *.  (?)] , 


(46) 


(47) 


which  is  Colby’s  result  when  his  parameter  6  is  evaluated  at  zero.20 

To  compare  Eqs.  (37)  and  (46),  Brillouin’s  choice  of  parameters  for  a  Lorentz  medium  are  used 
(u>o  =  4.0  X  1016s  *,  b2  =  20.0  X  1032s  2,  6  =  0.28  X  1016s-1),  and  the  frequency  uc  of  the  incident 
field  is  10  s  .  As  Fig.  8  illustrates,  at  an  observation  point  that  is  10-5  meters  to  the  left  of  the 
dispersive  medium’s  boundary,  E^et  (lower  curve)  and  Ery>h}  (upper  curve)  agree  extremely  well  for 
1  -  6  -  1-0013-  Bot]l  signals  are  not  observed  until  33.36 ’femtoseconds  (fs)  after  the  incident  field 
encounters  the  medium.  For  the  first  43.36  attoseconds  (as)  after  their  initial  observation,  they  are 
visually  indistinguishable.  At  that  point,  they  start  separating  rapidly  because  the  high-frequency 
approximation  is  no  longer  valid. 
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Fig.  8  -  Early-time  (lower  curve)  and  high-frequency  (upper  curve)  fields 
vs  normalized  time  9  =  ct/\z\  at  z  —  — 10-5  m  for  U>o  —  4.0  X  1016s_1, 
b 2  =  20.0  x  103V2,  and  6  =  0.28  x  K^V1. 


4.  EVALUATION  OF  TRANSMITTED  FIELD 

Since  a  closed-form  evaluation  of  the  inversion  integral  in  Eq.  (11(b))  is  unlikely,  one  must 
resort  either  to  numerical  evaluation  of  the  inversion  integral  or  to  an  asymptotic  approximation  of 
it.  In  this  report,  the  latter  approach  is  taken  for  observation  points  that  are  sufficiently  far  from 
the  dispersive  medium’s  boundary  to  ensure  the  validity  of  asymptotically  representing  the  fields. 
The  inversion  integral,  a  Generalized  Laplace  Integral  [21,22],  is  treated  with  existing  uniform 
asymptotic  methods.  Detailed  discussions  of  uniform  asymptotic  techniques  are  found  in  Refs. 
23-27. 


Details  of  the  methodology  for  obtaining  a  uniform  asymptotic  expansion  of  the  transmitted 
field  for  a  constant  (unity)  transmission  coefficient  are  found  in  Refs.  4  and  6-8.  The  purpose  of 
this  section  is  to  show  how  inclusion  of  a  frequency-dependent  transmission  coefficient  impacts  the 
transmitted  field. 

4.1  Inclusion  of  Transmission  Coefficient  in  Transmitted  Field 

In  the  notation  of  Oughstun,  Eq.  (11(b))  becomes 

/oo+ta 

A{u)T{u)eW“V  ^  (48) 

-oo+ia 

where  the  dimensionless  parameter  9  =  ctjz  can  be  viewed  as  normalized  time  and 

</>(u>,  9 )  =  iu[n(u>)  -  9 ].  (49) 

The  symbols  £  and  A  are  used  for  the  fields  when  they  depend  on  9  instead  of  time  t.  Although 
the  right  side  of  Eq.  (48)  is  similar  to  the  equations  given  in  Ref.  28,  it  differs  in  two  ways:  (1)  the 
missing  factor  of  27T,  which  comes  from  the  difference  in  the  definitions  of  the  Fourier  transform, 
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and  (2)  the  inclusion  of  the  transmission  coefficient  T{u>).  The  saddle  points  of  Eq.  (48)  are  the 
critical  points  of  the  generalized  phase  tf>,  which  are  the  roots  of 


n(u)  +  um'(w)  -5  =  0.  (50) 

When  5  <  1  (t  <  z/c),  one  can  show  that  £*(££)  =  0  by  closing  C\  with  a  semicircle  from 
above.  Thus  at  a  fixed  observation  point  z  within  the  dispersive  medium,  no  signal  is  observed  until 
the  leading  edge  of  the  incident  field  traverses  the  distance  z.  When  6  >  1,  it  is  well-known  that 
the  inversion  integral  can  be  evaluated  by  closing  C\  so  that  the  resulting  Jordan  contour  Cj(6) 
passes  through  the  saddle  points  via  the  subcontour  P(5)  =  Pn(0)  +  Ed(5)  [29,30]  and  so  that  the 
integrals  over  subcontours  other  than  C  and  P(0)  go  to  zero  via  some  limiting  process. 


To  evaluate  the  integral  representation  of  £*,  first  observe  that  the  Fourier  transform  M  of  an 
arbitrary  modulation  function  M  is  related  to  A  (Eq.  (10(a)))  by 


OJ  +  COc)  —  M((jJ  —  LOc 


>]• 


(51) 


The  difference  between  the  u  of  Oughstun  and  Sherman  and  M  arises  from  the  2jt  variance  in  the 
definitions  of  the  Fourier  transform;  that  is,  2 vM  =  u.  Substituting  Eq.  (51)  into  Eq.  (48)  yields 


—  1  T  „  roo+ia  -i 

£y(z’Q)  =  /  T{ui)M{u  +  uc)e‘^u’9^  dio  —  I  T(w)M{u  —  uc)e*^'e^  du  .  (52) 

z  L^-oo+*a  J-oo+ia  J 


Because  M  is  zero  for  t  <  0  and  is  real-valued  for  t  >  0,  M  satisfies  M(-oj)  =  [M(w*)]*.  Conse¬ 
quently,  <j>*{u*,9)  =  T*(u*)  =  T{-u),  and  M(-(u  -wc)*)  =  [M(w  -  wc)]*,  which  imply 

that  the  first  integral  in  Eq.  (52)  is  the  complex  conjugate  of  the  second  integral.  Hence, 


T  i  f°°+ia 

Sl(z,  6)  =  Re  —  T(u)u(u  -  ujc)e^^  dw 

LZ7r  J-oo+ia 


which  is  Eq.  (1.11)  of  Ref.  7  and  Eq.  (7.2.1)  of  Ref.  8  with  T(u>)  =  1. 


(53) 


To  this  point,  the  results  are  general.  When  M  is  the  window  function  defined  by  Eq.  (12), 

giorTo  _ 

i(")  =  — —  (54) 

and 

Si(z,9)  =  A(z,e,0)-A(z,$,To),  (55) 

where 


5t0  =  5 


cT0 


and  «4(z,5,T0)  =  Re 

27T 


(56) 


To  determine  the  transmitted  field,  one  need  only  derive  A(z,0,To)  and  take  the  limit  as  T0 
approaches  zero  to  obtain  .A(z,5,0).  For  -4(z,5,T0),  apply  Cauchy’s  residue  theorem  to  obtain  the 
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decomposition  of  Oughstun  and  Sherman, 


A(z,9,T0) 


e-tu/cT0 

/  ef0(w,«To)  fa 

27 r 

. 

IPo(6)  V  -U>c 

W w«. 

e-ta;cTo 

j  T(W)  e7<H“’eT0)  du 

27T 

Jp»(e )  w  ~  wc 

u(M 

2n 


Re 


—2 irie  tu>cT°  ^Residues  enclosed  by  Cj(9) 


As(z,0,Tq)  +  Ab(z,9,Tq)  +  Ac(z,9,To), 


(57) 


where  the  first,  second,  and  third  terms  following  the  first  equal  sign  correspond  to  As,  Ab,  and 
Ac,  respectively. 


Note  that  A(z ,  0,  To )  and  A(z ,  0, 0)  are  zero  for  6  <  1  and  0  <  cTo / z,  respectively,  and  that  each 
is  basically  the  integral  along  P(6 )  plus  contributions  from  poles  of  the  transmission  coefficient  T 
and  the  Fourier-transformed  modulation  function  M  that  lie  inside  Cj(6 ).  The  Generalized  Laplace 
Integrals  in  Eq.  (57)  are  evaluated  with  uniform  asymptotic  expansions  relative  to  the  secondary 
parameter  0  as  z  increases  without  bound  [23,25,26]. 

The  methodology  of  Oughstun  and  Sherman  [31],  which  combines  uniform  asymptotic  results 
from  five  sources  [27,32-35],  is  followed  to  get  a  uniform  asymptotic  representation  of  £*.  Defining 
Aj(z,0)  =  Aj{z,0, 0)  —  Aj(z,0,To)  for  j  in  {S',  B,C)  implies 

£*(*>  0)  -  As(z,  0)  +  Ab(z,  0 )  +  Ac(z,  9).  (58) 

The  contributions  As  (Sommerfeld’s  precursor  field),  Ab  (Brillouin’s  precursor  field),  and  Ac 
(main  signal)  depend  respectively  on  the  behavior  of  the  distant  saddle  points,  on  the  behavior  of 
the  near  saddle  points,  and  on  the  location  of  the  poles  of  T{u)M{ u  —  uc)  relative  to  path  P(0 ). 
Since  the  refractive  index  is  bounded  away  from  -1,  T  has  no  poles  and  is  bounded.  Hence,  the 
pole  contributions  come  solely  from  M,  whose  only  pole  is  located  at  u;c. 


4.2  Saddle  Points  of  <j> 


As  mentioned  earlier,  the  saddle  points  of  Eq.  (48)  are  the  roots  of  Eq.  (50).  After  substituting 
Eq.  (6)  into  Eq.  (50),  eliminating  the  square  root,  and  manipulating  the  resulting  equation,  one 
obtains  the  following  eighth-degree  polynomial  in  cj, 

0  =  [92  -  l]w8  +  [8tfi(02  -  l)]u>7  -  [4(6£2  +  ul){9 2  -  1)  +  b292]u6 

-  [26i({l2L>Z  +  362  +  16£2)(02  -  1)  +  2b2)]u5 

+  [(6o4  +  m2ul  +  2 b2ul  +  12 b2S2  +  16<54)(02  -  1)  +  62(c^02  -  12<52)  +  \6b2S2]u>4 
+  [4fa'((6u;o  +  Zb2ul  +  8 62ul  +  2b262){92  -  1)  -  b2(2ul  +  3<52)  +  b2( 2w\  +  3£2))]w3 

-  [wo(4a;o  +  3b2ul  +  2A82ul  +  12b262)(92  -  1) 

-  62(wq  +  2062u>1  +  9b2 S2)  +  8b262(2wl  +  b2)]u2 

-  [2 6iu>%  ((4 u>l  +  Zb2)(ojl(02  -  1)  -  62)  +  2b2(u>l  +  62))]a> 

+  [^o(wo  +  b2){u2(92  -  1)  -  62)] , 


(59) 
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that  must  be  solved  to  obtain  the  saddle  points.  It  turns  out  that  Eq.  (50)  has  exactly  four  saddle 
points.  Consequently,  the  squaring  operation  that  leads  to  Eq.  (59)  introduces  four  extraneous 
roots.  The  saddle  points  are  symmetric  with  respect  to  the  imaginary  axis  in  the  complex  u>  plane 
and  fall  naturally  into  two  categories,  those  that  are  near  the  origin  (u>m  and  uN2)  and  those  that 
are  far  from  the  origin  (u;£>i  and  U} D2). 

To  place  the  following  brief  discussion  in  context,  a  few  of  the  salient  features  of  the  dynamical 
behavior  (0  dependence)  of  the  saddle  points  [36]  are  first  summarized.  The  near  saddle  points  lie 
along  the  imaginary  axis  for  0  =  1,  with  uN1  being  the  upper  one.  As  0  increases,  they  approach 
each  other  on  the  imaginary  axis  until  they  coalesce  to  a  second-order  saddle  point  at  0\.  As  0 
continues  to  increase  beyond  01}  the  near  saddle  points  symmetrically  separate  with  increasing 
distance  between  them.  In  the  limit  as  0  -»•  oo,  they  approach  the  branch  points  u±.  On  the 
other  hand,  the  distant  saddle  points  always  lie  in  the  lower  half  of  the  complex  plane.  Moreover, 
they  are  located  at  ±oo  -  2 Si  when  0=1  and  move  towards  the  branch  cut  as  0  increases,  with 
lime-.oo  «di  =  u'+  and  lim^,*, u>D2  =  u>'_. 

On  comparing  Eq.  (59)  with  Eq.  (4.2.5)  in  Ref.  4  and  Eq.  (6.2.5)  in  Ref.  8,  one  observes 
discrepancies  m  the  coefficients  of  u1,  w2,  w3,  w4,  and  w5.  To  ascertain  the  impact  of  the  difference 
between  the  two  polynomials,  the  saddle  points  of  both  polynomials  are  calculated  for  Brillouin’s 
choice  of  Lorentz  parameters  (w0  =  4.0  x  1016/s,  b2  =  20.0  x  1032/s2,  S  =  0.28  x  1016/s).  These 
values,  computed  independently  with  Mathematica  and  MATLAB,  are  provided  in  Tables  3  and 
4.  Although  the  roots  are  computed  to  12  decimal  places,  the  tabulated  entries  are  rounded  to  3 
decimal  places  to  conserve  space.  The  values  of  0  at  which  the  roots  are  determined  axe  chosen 
so  that  the  computed  values  may  be  compared  with  values  that  are  visually  estimated  from  Fies 
6.15,  6.16,  and  6.17  of  Ref.  8. 
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Table  3  -  Distant  saddle  points  u>di,  roots  from  two 
distinct  eighth- degree  polynomials,  vs  normalized  time 
6.  Column  2  is  calculated  with  Eq.  (66).  Column  3  is 
calculated  with  Eq.  (6.2.5)  of  Oughstun  and  Sherman.a 


e 

U?D1  "  10  16 

(rad/s) 

udi  •  10"16 

(rad/s) 

1.000 

oo  -  0.560* 

oo  —  1.120t 

1.009 

34.245  -  0.541i 

34.264  -  1.062* 

1.010 

32.580  -  0.539* 

32.600  -  1.056* 

1.015 

26.971  -  0.530i 

26.993  -  1.029* 

1.020 

23.667  -  0.522* 

23.690-  1.005* 

1.025 

21.437  -  0.515i 

21.460  -  0.983* 

1.030 

19.807  -  0.508* 

19.830  -  0.962* 

1.040 

17.548  -  0.496* 

17.571  -  0.926* 

1.050 

16.032  -  0.486i 

16.054  -  0.895* 

1.100 

12.389  -  0.449* 

12.406  -  0.787* 

1.150 

10.849  -  0.427* 

10.863  -  0.719*' 

1.200 

9.960  -  0.411* 

9.971  -  0.672* 

1.250 

9.369  -  0.399* 

9.378  -  0.636* 

1.300 

8.940  -  0.389i 

8.948  -  0.608* 

1.400 

8.353  -  0.375*' 

8.359  -  0.565* 

1.500 

7.964  -  0.365* 

7.968  -  0.534* 

1.502735905 

7.955  -  0.364* 

7.959  -  0.533* 

1.510 

7.932  -  0.364* 

7.936  -  0.531* 

1.520 

7.901  -  0.363*' 

7.905  -  0.529* 

1.530 

7.871  -  0.362* 

7.874  -  0.526* 

1.540 

7.842  -  0.361* 

7.845  -  0.524* 

1.550 

7.814-  0.360*' 

7.817  -  0.521* 

1.600 

7.684  -  0.356* 

7.686  -  0.510*' 

1.700 

7.471  -  0.350* 

7.473 -0.490i  | 

1.800 

7.303  -  0.345* 

7.304  -  0.474*' 

1.900 

7.166  -  0.340* 

7.167  -  0.460* 

2.000 

7.053  -  0.336*' 

7.053  -  0.448* 

2.500 

6.690  -  0.322* 

6.688  -  0.405* 

3.000 

6.493  -  0.313*' 

6.491  -  0.378* 

4.000 

6.290  -  0.302*' 

6.287  -  0.346*' 

5.000 

6.190  -  0.295* 

6.187  -  0.327* 

10.000 

6.046  -  0.285* 

6.045  -  0.294* 

20.000 

6.007  -  0.281* 

6.007  -  0.284*' 

100.000 

5.994  -  0.280* 

5.994  -  0.280* 

1000.000 

5.993  -  0.280* 

5.994  -  0.280* 

oo 

5.993  -  0.280* 

5.993  -  0.280*' 

aK.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propa¬ 
gation  in  Causal  Dielectrics  (Springer- Verlag,  Berlin,  1994),  Ch.  6, 

p.  208. 
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Table  4  Near  saddle  points  {u>m^N2},  roots  from  two  distinct  eighth-degree  polynomials,  versus 
normalized  time  6.  Columns  2  and  4  are  calculated  with  Eq.  (66).  Columns  3  and  5  are  calculated 
with  Eq.  (6.2.5)  of  Oughstun  and  Sherman* 


6 

U)n\  *  10  16 
(rad/s) 

&N1  *  10  16 
(rad/s) 

<^N2  ’  10  16 

(rad/s) 

Wj\T2  •  10  16 
(rad/s) 

1.000 

4.474  i 

5.080  ? 

-4.461? 

-3.924? 

1.009 

4.297  ? 

4.854  ?' 

-4.322? 

-3.814? 

1.010 

4.278  i 

4.831  ?' 

-4.307? 

-3.802? 

1.015 

4.189  ? 

4.720  ?' 

-4.236? 

-3.745? 

1.020 

4.105  ?' 

4.616  ? 

-4.168? 

-3.690? 

1.025 

4.024  i 

4.518  i 

-4.102? 

-3.636? 

1.030 

3.948  ? 

4.426  ?' 

-4.039? 

-3.584? 

1.040 

3.805  ? 

4.255  ? 

-3.920? 

-3.485? 

1.050 

3.673  i 

4.100  ?' 

-3.809? 

-3.391? 

1.100 

3.127  ?' 

3.474  ? 

-3.335? 

-2.981? 

1.150 

2.695  i 

2.995  ?' 

-2.948? 

-2.635? 

1.200 

2.326  i 

2.592  ? 

-2.610? 

-2.326? 

1.250 

1.992  i 

2.234  ? 

-2.300? 

-2.038? 

1.300 

1.677  i 

1.900  ?' 

-2.003? 

-1.758? 

1.400 

1.040  i 

1.230  ? 

-1.393? 

-1.171? 

1.500 

0.000  i 

0.000  -  0.000  ?' 

-0.373? 

-0.000? 

1.503b 

O 

o 

o 

o 

1 

o 

►-» 

00 

0.186  -  0.001? 

© 

o 

o 

o 

1 

© 
J— » 
00 

0.186  -  0.001? 

1.510 

0.302  -  0.187? 

6.355  -  0.003? 

0.302  +  0.187? 

0.355  -f*  0.003t 

1.520 

0.463  -  0.188? 

0.499  -  0.005? 

0.463  +  0.188? 

0.499  +  0.005? 

1.530 

0.579  -  0.189? 

0.608  -  0.007? 

0.579  +  0.189? 

0.608  +  0.007? 

1.540 

0.673  -  0.190? 

0.698  -  0.010? 

0.673  +  0.190? 

0.698  +  0.010? 

1.550 

ro.754  -  0.191? 

0.776  -  0.012? 

0.754  +  0.191? 

0.776  +  0.012? 

1.600 

r  1.054  -  0.194?' 

1.069-0.023? 

1.054  +  0.194? 

1.069  +  0.023? 

1.700 

1.432  -  0.200? 

1.442  -  0.041? 

1.432  +  0.200?' 

1.442  +  0.041? 

1.800 

1.685  -  0.205? 

1.693  -  0.056? 

1.685  +  0.205? 

1.693  +  0.056? 

1.900 

1.875  -  0.209? 

1.882  -  0.068? 

1.875  +  0.209? 

1.882  +  0.068? 

2.000 

2.026  -  0.213? 

2.031  -  0.079? 

2.026  +  0.213? 

2.031  +  0.079? 

2.500 

2.490  -  0.225? 

2.493  -  0.116? 

2.490  +  0.225? 

2.493  +  0.116? 

3.000 

2.742  -  0.233? 

2.745  -  0.139? 

2.742  +  0.233? 

2.745  +  0.139? 

4.000 

3.025  -  0.242? 

3.027  -  0.166? 

3.025  +  0.242? 

3.027  +  0.166? 

5.000 

3.186  -  0.248? 

CO 

U* 

oo 

1 

p 
► — k 
oo 

CO 

N. 

3.186  +  0.248? 

3.187  +  0.183? 

10.000 

3.510  -  0.260? 

3.511  -  0.219? 

3.510  +  0.260? 

3.511  +  0.219? 

20.000 

3.694  -  0.267?' 

3.695  -  0.241? 

3.694  +  0.267? 

3.695  +  0.241? 

100.000 

3.890  -  0.275? 

3.890  -  0.266? 

3.890  +  0.275? 

3.890  +  0.266? 

1000.000 

3.969  -  0.279?' 

3.969  -  0.277? 

3.969  +  0.279? 

3.969  +  0.277? 

00 

3.990  -  0.280? 

3.990  -  0.280? 

3.990  +  0.280? 

3.990  +  0.280? 

aK-E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994),  Ch.  6,  p.  208. 

1.502735905  is  rounded  to  conserve  space. 
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The  values  calculated  with  Eq.  (49)  (column  2  of  Table  3,  columns  2  and  4  of  Table  4)  do 
not  match  those  calculated  with  Eq.  (6.2.5)  of  Ref.  8  (column  3  of  Table  3,  columns  3  and  5  of 
Table  4).  Specifically,  the  imaginary  parts  of  the  distant  saddle  points  differ  considerably.  Until  6 
exceeds  2,  the  imaginary  parts  are  separated  by  at  least  0.1.  In  fact,  the  greatest  disparity  occurs 
near  0  =  1,  where  apparently  the  roots  of  Oughstun ’s  polynomial  approach  oo  -  1.1200i  =  oo  —  46i. 
This  value  is  estimated  by  determining  the  roots  for  1  <  6  <  1.009.  Clearly,  this  disagrees  with 
the  theoretically  predicted  value  of  oo  -  2 Si]  whereas  from  Eq.  (59),  the  limiting  value  of  upx 
as  0  — ►  1+  is  in  agreement  with  the  theoretical  prediction.  Similarly,  for  the  near  saddle  points 
(up? i  and  0^2),  the  difference  between  the  roots  is  largest  at  6  =  1  and  gradually  diminishes 
as  0  increases.  As  Oughstun  and  Sherman  note,  uoi  -*  =  (5.9935  -  0.280(H)  x  1016  and 

uni  —►«*>+  =  (3.9902  -  0.2800i)  x  1016  as  6  approaches  infinity.  In  addition,  the  critical  value 
where  the  near  saddle  points  change  from  pure  imaginary  to  being  symmetrically  located  on  either 
side  of  the  imaginary  axis,  is  1.502735905  to  eight  decimal  places. 

Although  the  numerical  evaluations  of  the  roots  via  Eq.  (59)  differ  from  roots  that  are  computed 
with  the  eighth-degree  polynomial  of  Oughstun  and  Sherman,  they  do  match  the  figures  on  pages 
236-238  of  Ref.  8.  Consequently,  the  authors  suspect  that  in  calculating  the  “exact”  roots  for  their 
figures,  Oughstun  and  Sherman  applied  Newton’s  method  to  an  equation  other  than  the  polynomial 
given  by  Eq.  (6.2.5). 

4.3  Contributions  from  Distant  Saddle  Points 


In  this  section,  a  uniform  asymptotic  approximation  of  As{z,0)  is  obtained  by  following  the 
derivations  in  Ref.  7  and  Section  7.2  of  Ref.  8.  From  Theorem  5.4  of  Ref.  8,  a  restatement  of  a 
theorem  in  Ref.  27,  a  uniform  asymptotic  expansion  of  As(z,0^To)  to  first  order  is 


As(z,e,T0)  = 


o-tu>cT0 


U(eT b)Re 


/  o)duj  (60a) 

JPdWu~uc 

e-,z(3(6To)  [2 a(0ro)e~,^]1/  7o(0To )Jv(za(6To)) 

-t-2a(0To)e-,^7i(0To)Jl/+i(zQ:(0ro))j| ,  (60b) 


where  Jv  and  Ju+\  are  Bessel  functions, 
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<*(*Xo  )  =  Yc  [^DX  i  @T0  )  -  ^(W£)2 ,  Ox )] , 
0(@To  )  =  2^  [^(^£>1 ,  #To  )  +  4>(uD2  ,  Oxo  )] » 


7o(9t-) = ns(“o,(9T')  -  “e)  [4d  ,+4-<m  [&((£:UJ 5 

+  .(«„<*)  -  «.)  [3^]%^))  ’}. 

7,(#R)  =  {*“ "  “J  [4s]  1+V“»  [i^^k) 

-  «•>  -  «■>  [s^S)k)  1  ’  )• 


The  subscript  u>u  of  <j>  denotes  the  second  partial  derivative  of  <f>  with  respect  to  u>,  and  v  is  a 
real  number  such  that  T(u>)/(u>  —  u>c)  =  u)~^1+^g(u>)  and  g  has  a  convergent  Laurent  series  about 
w  =  0.  For  the  window  function  of  Eq.  (12),  v  =  0.  The  function  v  corresponds  to  u  in  Eqs. 
(7.2.13-15)  of  Ref.  8.  When  one  sets  T(ojDj(0To))  =  1  and  identifies  a  and  0  with  a/c  and  0/c  of 
Ref.  8,  Eqs.  (60)  and  (61)  reduce  to  Eqs.  (7.2.11-15)  of  Ref.  8.  The  reason  for  this  slight  difference 
in  nomenclature  is  that  Oughstun  and  Sherman  normalize  the  a  and  0  of  Theorem  5.4  [8]  with 
division  by  c.  Regardless,  the  final  results  are  equivalent  when  the  transmission  coefficient  is  unity. 

The  approach  taken  in  Section  7.2  of  Ref.  8  approximates  the  expressions  in  Eqs.  (61)  directly; 
however,  Theorem  5.4  actually  provides  some  asymptotic  approximations  for  a,  0,  and  the  radical 
in  Eqs.  (61).  Before  proceeding,  the  authors  believe  that  Eqs.  (5.4.10-12)  of  Ref.  8  and  Eqs.  (3.4.10- 
12)  of  Ref.  4  are  not  quite  correct  as  a  consequence  of  errors  in  some  works  of  Handelsman  and 
Bleistein  [37].  These  errors  have  been  corrected  on  page  390  of  Ref.  25  and  result  in  the  following 
modifications  to  expressions  that  occur  in  pertinent  works  of  Oughstun  and  Sherman, 


a(0)  =  2[-ai(0)(e  -  1)]*  +  o{[6  -  l)i); 

m  =  ao(0)  +  0{(0-l)), 


(62a) 

(62b) 


4  az(0) 
>(2)(u,±,0) 


:4M0)||l  +  o((0-l)*)  . 


As  Eqs.  (61)  indicate,  one  derives  an  approximation  to  «4s(z,  0,To)  by  developing  accurate 
estimates  of  the  distant  saddle  points  and  generalized  phase.  According  to  the  so-called  second 
approximation  to  the  distant  saddle  points  [38,39], 

WDi(0T„)  =  £(07o)  -  w[l  +  (63a) 

W£>2(0T„)  =  -f(0To)  -  ^[1  +  f?(0ro)],  (63b) 

=  “(I-*)- 2(^4? 


(63c) 
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where 


b2$T  1 1 

£{@To)  =  |«0  ^  +  ^2  _  ij  > 

(64a) 

ri{eTo)  =  i 

(64b) 

Substituting  Eqs.  (63)  and  (64)  into  Eqs.  (60)  and  (61)  leads  to  more  general  versions  of 
7o,  71,  and  As  in  Eqs.  (7.2.13-15)  of  Ref.  8.  In  hindsight,  one  can  obtain  these  generalizations  by 
replacing  u{usp+  — wc)  and  u{usp-  -uc)  of  Ref.  8  with  v(u£>i  —vc)T(<jJdi  )  and  v{udi  — u>c)T(u>/)2)> 
respectively.  Let  Tr,  Tj,  tir,  and  n;  be  the  real  and  imaginary  parts  of  the  transmission  coefficient 
(T  =  Tr  +  iTi)  and  the  refractive  index  (n  =  hr  +  in/);  then  by  Eq.  (20a), 


T  (i.A  -  2(1  +  ngH) 

'  (l  +  nH(u;))2  +nj(u) 

T  (t.A  -  ~2n/H 

^  ^  (1  +  hk(w))2  +  n2j(u) 


«n(wDi,2(^T0))  =  1  - 


62r2(^o)[l  +  ^r2(gTo)[l-7?2(gTo)]] 

2[l  +  ^-2(0To)[1  +  ^To)]2]  [l  +  ^-2(^T„)[1  -  V(0tb)]2}  ’ 


W/(^£)l,2(^To)) 


_ 3{8t0)v(6t0) _ 

2[l  +^-2(^t„)[1  +  ^T„)]2]  [l  +  ^-2(0to)[1  -  v(0To)}2'  ' 


(65a) 

(65b) 

(65c) 

(65d) 


The  dependencies  of  the  right  sides  of  Eqs.  (65(a)  and  (b))  on  w  are  not  shown  explicitly  but  are 
specified  by  Eqs.  (19).  In  comparing  Eqs.  (65(a)  and  (b))  and  Eq.  (20(b)),  note  that  the  imaginary 
parts  of  T( u)  and  i?(u>)  are  equal  and  that  T(u>)  —  R(oj)  =  Re[T(w)  —  R(u>)]  =  1.  Moreover,  the 
approximations  for  tir  at  u>m(QT0)  and  udi(Qt0)  are  equal;  however  the  minus  and  plus  signs  in 
Eq.  (65(d))  correspond  to  n/  evaluated  at  u>di(^t0)  and  wD2(^t0)>  respectively. 
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Finally,  by  substituting  Eqs.  (63)  and  (65)  into  Eqs.  (60)  and  (61),  one  obtains 


As(z,0,To)  =  U(eTo)^l^To- 


1  + 


b2 


13 


2mo)+^[i-^To)pj 


x  e 


x 


+ 


+ 


+ 


X  J0 


2H(^JD1(^To))[f^e(^T0)  -  wc][l  -  rj(9To)]  +  S^(0To)[l  +  J?(0To)]l 
[^T0)-«cP+^[l  +  7?(0To)P 
Ti(u>d i(9To)) [^To)[^(^To)  -  Wc]  -  §<52[1  -  7p(0To)]] 
[^To)-a;cp  +  ^[l  +  7?(tfTo)P 

7fi(a;I5i(eTo))[|^[^To)  +  Wc][l  -  t)(6t0 )]  +  ££(0to)[1  +  v{^Ta )]] 

[^To)+^P  +  ^[1  +  ^To)]2 

7’/(wr>l(^To))[^To)[^(^To)  +  wc]  -  §<52[l  -  7p(0To)]j  \ 

[?(^To)  +  WcP  +  <?2[1  +  n(^To)]2  J 

- 1+g  1  ]) 

1  2^To)+^[i-^to)pJ; 

^7i?(wDx(^To))[^To)[^To)  -  Wc]  -  §(52[1  -  J?2(0To)]| 

K(*  To)-WcP+^[l  +  7?(0To)p 

I/(wjDi(Oro))[f^(0T„)  -  u>c][l  -  v(6To)]  +  ^(0To)[l  +  j/(0Xo)]] 

[^To)-^P+^[l  +  7/(0To)P 

_  TR(uD1(6To)) [^To)[f(^o)  +  -  f<52[l  -  7?2(gTo)3] 

[C(^T0  )  +  Wcp  +  £2[1  +  7?(0To)P 

r/(wDl(^T0))[|^(^To)  +  wc][l  -  V(8t0  )]  +  ^(*T„)[1  +  7?(0TO)]|  > 
[^To)  +  a;CP+^2[l  +  ^To)P  i 


V 


+ 


X  JA 


1  + 


2^2(0to)  +  ^[1-^t„)PJ 


))' 


(66) 


In  the  limit  a,s  TR  ->  1  and  Tj  -»  0,  As(z,6,  T0)  reduces  to  Eq.  (7.2.24)  of  Ref.  8.  The  effect  of  T 
on  ^4 s  is  negligible  because  T  =  1  for  the  extremely  high  frequencies  associated  with  this  precursor. 


4.4  Example  of  Transmitted  Field 


The  contributions  As(z,0)  to  the  transmitted  field  from  the  distant  saddle  points  {u)D1,ojD2} 
are  evaluated  and  plotted  for  the  observation  point  at  z  -  10-5m  in  the  Lorentz  medium,  for 
Brillouin’s  values  of  the  physical  constants  (u;0  =  4.0  x  1016s-1,  b2  =  20.0  x  1032s-2,  6  =  0.28  x 
1016s  2),  and  for  an  incident  plane  wave  that  consists  of  a  single  cycle  of  the  sine  at  angular  carrier 
frequency  we  =  1016s  2.  With  this  choice  of  parameters,  the  medium  is  highly  absorptive  about 
u>0)  and  uc  (To  =  2w/uc)  is  below  the  medium’s  absorption  band,  ((wq  —  S2)1^2 , (a;2  —  tf2)1/2). 
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The  results  for  As  are  then  compared  to  the  contributions  given  by  Eq.  (7.2.24)  of  Ref.  8,  where 
Oughstun  and  Sherman’s  representation  of  the  Sommerfeld  precursor  is  designated  Aso{z,9). 

Graphs  of  As(z ,  9 , 0)  and  As(z,  9,  To)  versus  normalized  time  9  are  shown  in  Figs.  9(a)  and  (b), 
respectively.  As  one  expects  from  Eqs.  (49),  (56),  and  (60),  As(z,9, To)  is  As(z,9, 0)  delayed  by 
cTo/z  =  0.0188365.  On  inspection,  the  plot  of  As(z,9, 0)  looks  like  Fig.  6  of  Ref.  8.  By  substituting 
To  =  0  in  Eq.  (66)  and  by  manipulating  Eqs.  (63-66),  one  can  argue  that  the  envelope  of  As(z,  9, 0) 
in  Fig.  9(a)  behaves  like 

±(9  -  (67a) 

Note  that  the  exponent  has  linear  and  quadratic  terms  in  (9  —  1).  For  9  near  1,  the  exponent  is 
dominated  by  the  linear  term,  and  Eq.  (67a)  becomes 

±(0- l)e-(2z5/c)(9-1),  (67b) 

which  is  proportional  to  the  Gaussian  probability  density  function  with  parameters  2  and  c/{2z6) 
[40].  On  the  other  hand,  when  9  is  large,  the  envelope  behaves  like  the  Rayleigh  probability  density 

±(9-  i)e-0«/(2c))(«-i)2.  (67c) 


Recall  that  As(z,  9)  =  .4s(z,0,O)  -  As(z,9,To).  When  this  difference  is  plotted,  one  observes 
two  slightly  overlapped  signal  packets  (Fig.  9(c)),  which  correspond  to  the  intervals  1  <  9  <  cT^/z 
and  9  >  cTo/z.  In  the  first  packet  (1  <  9  <  cTo/z),  As(z,9 )  =  As(z,9,0)  since  As(z,9,To)  =  0. 
Moreover,  at  9  =  cTo/z,  the  magnitude  of  As(z,9)  has  decreased  roughly  to  one  tenth  of  its 
maximum  value.  When  9  >  cTo/z,  the  contribution  from  As(z,9, 0)  decreases  as  9  increases. 
However,  until  .As(z,0,O)  becomes  negligible  near  9  =  1.04,  the  trailing  edge  of  the  first  packet 
interferes  with  the  leading  edge  of  the  second  packet,  thereby,  producing  the  uneven  behavior  of 
the  second  packet’s  envelope.  On  comparing  Fig.  9(c)  to  the  first  part  of  the  bottom  frame  of  Fig. 
15  in  Ref.  7,  the  curves  appear  to  be  very  similar.  To  ascertain  the  degree  of  similarity  between 
As  and  Ago,  the  difference  As{z,9)-Aso(z,9)  is  plotted  in  Fig.  9(d).  As  the  figure  indicates,  the 
absolute  difference  is  always  less  than  0.000004.  Thus  the  graph  of  As  is  visually  indistinguishable 
from  Aso •  The  reason  for  this  lies  in  the  behavior  of  T(u>£n  (9)).  According  to  Table  A1  of  the 
Appendix, 

1  <  TK(w0i(0))  <  1.02  and  0  <  T/(ua(0))  <  0.00062,  (68) 

for  1  <  9  <  1.05,  which  implies  that  Tr  =  1  and  T/  =  0  at  values  of  9  for  which  the  Sommerfeld 
precursor  is  significant.  On  the  other  hand,  the  decay  of  the  exponential  factor  completely  domi¬ 
nates  the  contribution  from  the  variations  in  T(ujx(9))  when  9  exceeds  1.05  since  T  is  bounded. 
Hence  choosing  T(udl(9))  A  1  in  Eq.  (66)  has  a  negligible  effect,  and  the  small  excursions  of  the 
transmission  coefficient  from  unity  at  optical  frequencies  for  Brillouin’s  constants  have  little  impact 
on  As-  This  analysis  merely  validates  the  tacit  assumption  that  T  =  1  at  optical  frequencies  for  the 
Sommerfeld  precursor.  Of  course,  one  expects  the  impact  of  T(u>)  on  As  to  be  different,  perhaps 
more  significant,  over  other  frequency  regimes  like  millimeter  waves  or  microwaves. 
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Fig.  9  -  Sommerfeld  precursor  and  associated  functions  versus  normalized  time  6  at  z  =  -10~5  m  in  the  Lorentz 
medium  with  u0  =  4.0  X  lO16®*1,  b 2  =  20.0  X  1032s"2,  and  S  =  0.28  X  lO'V1:  (a)  AS(z,0, 0),  (b) 
•As(z’  To),  (c)  As(z,  0)  =  As(z,  0, 0)  -  As{z,  0,  To),  and  (d)  As(z,  0)  -  Aso(z,  0).  The  incident  pulsed 
field  is  a  single-cycle  sine  at  U>c  =  lO^s"1  with  T0  =  2tt/uc,  and  AS0  is  AS  evaluated  at  Tr(oj>)  =  1  and 
Tj(u>)  =  0.  v  7 


As  stated  in  the  introduction,  the  lengthy  endeavor  of  deriving  As  and  Ac  for  all  eventualities 
is  not  addressed.  However,  some  comments  are  made  regarding  the  impact  of  the  transmission 
coefficient  on  AB  and  Ac.  By  Eqs.  (57)  and  (58),  the  inclusion  of  T(u)  affects  the  amplitudes  of 
Ab  and  Ac  and  not  their  durations.  By  extrapolating  from  Figs.  8  and  10  of  Ref.  7,  one  deduces 
that  the  supports/peak  magnitudes  [41]  of  Ab  and  Ac  are  approximately  1.43  <  0  <  1.62/0.12  and 
1.44  <  0  <  1.70/0.005,  respectively.  Further,  by  Fig.  A2  and  some  calculations,  As  is  negligible 
for  0  >  1.08.  Consequently,  the  transmitted  field  has  a  significant  nonzero  value  for  0  in  [1,1.08], 
[1.43, 1.62],  or  [1.44, 1.70].  Although  further  observations  on  Ab  and  Ac  can  be  made  only  after 
detailed  study,  one  more  conclusion  regarding  Ac  can  be  drawn.  Since  T(uic)  =  0.7914  -  0.0051i 
is  a  multiplicative  factor  in  the  residue 


T{Uc)eWcW“‘’eT o)  (69) 

of  Ac  and  since  T(wc)  differs  substantially  from  unity,  the  inclusion  of  T(uc)  in  evaluating  Ac  will 
result  in  a  21%  decrease  in  the  main  signal’s  amplitude  as  computed  by  Ref.  7. 
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To  compare  As  to  Ey ,  consider  Fig.  6(d),  which  depicts  Ey  at  an  equal  distance  from  the 
dispersive  medium.  The  reflected  field  is  essentially  nonzero  for  0  in  [1,1.03]  and  has  a  peak  value 
of  0.21016  at  6  =  1.0142,  which  is  significant  relative  to  the  derived  and  expected  components 
of  Ey,  According  to  Table  2,  this  corresponds  nicely  to  the  value  of  |Re[iZ(cc?c)]  |.  Furthermore, 
|Re[R(u;c)]|  is  consistent  with  the  predicted  21%  percent  decrement  of  Ac •  Intuitively,  one  expects 
this  because  T(u;c)  -  R(u>c)  =  1.  The  comparison  between  the  expected  strengths  and  durations  of 
Ey  and  Ey  in  this  example  provides  a  compelling  argument  for  including  the  frequency  dependence 
of  the  reflection  and  transmission  coefficients  in  formulating  this  propagation  problem,  even  for 
fields  at  optical  frequencies. 

5.  SUMMARY  AND  CONCLUSIONS 

This  report  generalizes  the  formulations  of  Sommerfeld,  Brillouin,  and  Oughstun  and  Sherman 
for  one-dimensional  propagation  of  pulsed  fields  at  optical  frequencies  in  a  medium  for  which  the 
Lorentz  model  is  valid  and  for  which  the  Lorentz  absorption  frequency  6  is  less  than  its  resonance 
frequency  wo-  The  generalizations  are:  (1)  providing  a  formulation  that  is  applicable  to  frequencies 
for  which  the  Lorentz  model  of  the  medium  is  appropriate;  (2)  accomodating  a  pulsed  source  outside 
the  Lorentz  medium  instead  of  specifying  the  field  on  the  medium’s  boundary;  (3)  accounting  for 
reflection  from  the  medium;  and  (4)  including  a  frequency-dependent  transmission  coefficient  in  the 
transmitted  field.  Specifically,  results  are  obtained  for  the  case  when  the  source  is  a  rectangularly 
modulated  sine  function  at  angular  carrier  frequency  ljc  and  when  the  duration  To  of  the  source 
is  a  positive  integer  multiple  of  2i r/a>c.  Moreover,  results  on  the  reflected  and  transmitted  fields 
are  obtained  for  Brillouin’s  choice  of  Lorentz  parameters  (u>o,6,b).  This  example  demonstrates 
the  necessity  of  including  the  transmission  and  reflection  coefficients  for  this  class  of  propagation 
problems. 

An  exact  expression  for  the  time-domain  reflected  field,  that  is  amenable  to  physical  inter¬ 
pretation  and  numerical  evaluation,  is  derived.  When  time  t  =  0  is  designated  as  the  instant 
the  incident  field  reaches  the  medium’s  boundary,  the  reflected  field  is  first  observed  at  a  point  a 
distance  \z\  from  this  boundary  at  time  \z\/c.  The  field  consists  of  two  components:  (1)  a  sine 
function,  |J?(u>c)|  sm(uct  —  arg[R(u>c)]),  which  comes  from  the  pole  contributions  to  the  inversion 
integral  and  has  the  same  duration  To  as  the  incident  field  ( R  is  the  reflection  coefficient);  and  (2) 
two  exponentially  damped  sinusoidal  terms,  which  come  from  the  branch-cut  contributions,  with 
time- dependent  phases  that  broaden  the  incident  pulse  and  cause  ringing  in  the  trailing  edge.  The 
behavior  of  the  reflected  field  changes  with  the  value  of  uc  relative  to  u>o-  As  uc  increases,  the  field 
evolves  from  a  reduced,  negative  replica  of  the  incident  field  for  u>c  <  o>o/10  to  a  much  broader  field 
with  an  increasingly  distorted  leading  edge  and  a  damped  trailing  edge  having  a  greater  number 
of  zero  crossings  at  the  higher  carrier  frequencies.  For  uc  <  3u?o,  the  reflected  field’s  amplitude  is 
significant,  even  for  optical  frequencies.  In  addition,  early-time  and  high-frequency  fields  are  de¬ 
veloped,  and  they  are  shown  to  be  in  excellent  agreement  over  a  short  time  for  a  particular  optical 
example.  The  high-frequency  field  extends  Colby’s  result  by  including  the  attenuation  factor 
which  also  adds  higher-order  temporal  terms.  The  inclusion  of  these  terms  improves  the  accuracy 
of  Colby’s  result. 

By  following  the  methodology  of  Oughstun  and  Sherman,  the  time-domain  transmitted  fields 
are  expressed  as  the  sum  of  three  components:  the  Sommerfeld  precursor,  the  Brillouin  precursor, 
and  the  main  signal.  Rather  than  derive  uniform  asymptotic  representations  of  all  three  compo¬ 
nents  of  the  transmitted  field,  a  uniform  asymptotic  approximation  of  the  Sommerfeld  precursor  is 
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developed  to  illustrate  the  impact  of  the  transmission  coefficient  on  the  field.  Moreover,  the  Som- 

merfeld  precursor  is  shown  to  reduce  to  the  results  of  Oughstun  and  Sherman  when  the  transmission 
coefficient  is  unity. 

The  derived  theoretical  results  are  made  more  concrete  by  applying  them  for  Brillouin’s  choice 
of  the  Lorentz  medium’s  parameters  and  for  Oughstun ’s  choice  of  carrier  frequency.  Specifically,  the 
reflected  field  and  Sommerfeld  precursor  are  obtained  at  distance  \z\  from  the  medium’s  boundary. 
When  the  incident  pulsed  field  is  a  single  cycle  of  the  sine  function  (T0  =  2 k/uc),  neither  field  is 
observed  until  t  exceeds  |z|c-1.  The  Sommerfeld  precursor  has  a  significant  value  for  |z|c-1  <  t  < 
|z|c_1  +  47b  and  consists  of  two  slightly  overlapped  wave  packets.  The  envelope  of  each  packet  has 
a  leading  edge  that  resembles  the  Gamma  probability  density  function,  has  a  trailing  edge  similar 
to  the  Rayleigh  density  function,  and  bounds  a  highly  oscillatory  field.  The  transmission  coefficient 
has  a  negligible  effect  on  the  Sommerfeld  precursor  because  it  consists  of  extremely  high  frequencies; 
however,  a  heuristic  analysis  implies  that  its  impact  on  the  main  signal  is  significant.  In  contrast’ 
the  reflected  field  is  a  distorted  version  of  the  negative  of  the  incident  field  and  has  an  amplitude 
that  is  21%  of  the  incident  field’s  amplitude  and  is  roughly  250  times  the  peak  of  the  Sommerfeld 
precursor’s  envelope.  Further,  the  reflected  field’s  support  essentially  is  |z|c-1  <  t  <  \z\c~1  -f  1.5To, 
which  is  37.5%  of  the  precursor’s  support.  Given  these  facts,  one  is  drawn  to  the  conclusion: 
accurate  analyses  of  problems  involving  propagation  of  a  pulsed  field  that  is  externally  incident  on 
a  Lorentz  medium  require  frequency- dependent  reflection  and  transmission  coefficients. 

6.  ACKNOWLEDGMENT 

The  authors  thank  F.C.  Lin  for  providing  some  computations  with  MATLAB.  This  work  was 
sponsored  by  the  Naval  Research  Laboratory  and  the  Office  of  Naval  Research. 

REFERENCES 

1.  A.  Sommerfeld,  Ann.  Phys.  (Leipzig)  44(10),  177-202  (1914). 

2.  L.  Brillouin,  Ann.  Phys.  (Leipzig)  44(10),  203-240  (1914). 

3.  L.  Brillouin,  Wave  Propagation  and  Group  Velocity  (Academic  Press,  New  York,  1960). 

4.  K.E.  Oughstun,  “Propagation  of  Optical  Pulses  in  Dispersive  Media,”  Doctoral  dissertation, 
The  University  of  Rochester,  Rochester,  N.Y.,  1978. 

5.  K.E.  Oughstun  and  G.C.  Sherman,  J.  Opt.  Soc.  Am.  B  5(4),  817-849  (1988). 

6.  K.E.  Oughstun  and  G.C.  Sherman,  J.  Opt.  Soc.  Am.  A  6(9),  1394-1420  (1989). 

7.  K.E.  Oughstun  and  G.C.  Sherman,  Phys.  Rev.  A  41(11),  6090-6113  (1990). 

8.  K.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994). 

9.  A  radar  is  defined  to  be  ultrawideband  (UWB)  when  its  fractional  bandwidth  exceeds  25%, 
according  to  the  executive  summary  of  the  OSD/DARPA  study,  Assessment  of  Ultra-Wideband 
(UWB)  Technology ,  in  the  IEEE  Aer.  and  El.  Mag.,  November,  1990,  pp.  45-49. 

10.  B.  Noel,  ed.,  Ultra-Wideband  Radar:  Proceedings  of  the  First  Los  Alamos  Symposium  (CRC, 
Boca  Raton,  Fla.,  1991). 


34 


Mokole  and  Samaddar 


11.  H.L.  Bertoni,  L.  Caiin,  and  L.B.  Felsen,  eds.,  Ultra-Wideband,  Short-Pulse  Electromagnetics 
(Plenum,  New  York,  1993). 

12.  L.Yu.  Astanin,  A.A.  Kostylev,  Yu.S.  Zinoviev,  and  A.Ya.  Pasmurov,  Radar  Target  Character¬ 
istics:  Measurements  and  Applications  (CRC,  Boca  Raton,  Fla.,  1994). 

13.  J.D.  Taylor,  ed.,  Introduction  to  Ultra-Wideband  Radar  Systems  (CRC,  Boca  Raton,  Fla., 
1995). 

14.  L.  Carin  and  L.B.  Felsen,  eds.,  Ultra-Wideband,  Short- Pulse  Electromagnetics  2 (Plenum,  N.Y., 
1995). 

15.  W.  Colby,  Phys.  Rev.  V(3),  253-265  (1915). 

16.  A.  Sommerfeld,  Ann.  Phys.  (Leipzig)  44(10),  190  (1914). 

17.  L.  Brillouin,  Wave  Propagation  and  Group  Velocity  (Academic  Press,  N.Y.,  1960),  Ch.  Ill,  pp. 
56-57;  and  K.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal 
Dielectrics  (Springer- Verlag,  Berlin,  1994),  Ch.  6,  p.  234. 

18.  The  fractional  bandwidth  Bp  of  the  single-cycle  sine  function  is  given  by  the  formula,  Bp  = 
[2((jJu  —  wl)/(u>u  +  up )]  x  100%,  where  u>u  and  ul  are  the  upper  and  lower  frequencies  of  the 
passband  of  the  power  spectral  density  of  the  single-cycle  sine,  4u>2  sin2(7rw/u>c)/(«2  —  a;2)2. 
If  one  uses  the  half-power  criterion  to  determine  the  passband,  then  lou  =  1.309wc  and  up  = 
0.411u>c.  Hence  Bf  =  104%,  which  is  ultrawideband  according  to  Ref.  9. 

19.  L.L.  Pennisi,  Elements  of  Complex  Variables  (Holt,  Rinehart,  and  Winston,  New  York,  1963), 
p.  185. 

20.  W.  Colby,  Phys.  Rev.  V(3),  264  (1915). 

21.  A.  Erdelyi,  Asymptotic  Expansions  (Dover,  New  York,  1956),  Ch.  II,  pp.  29  and  36. 

22.  F.W.J.  Olver,  Asymptotics  and  Special  Functions  (Academic  Press,  New  York,  1974),  Ch.  3, 

p.  80. 

23.  E.T.  Copson,  Asymptotic  Expansions  (Cambridge,  London,  1965),  Ch.  10,  pp.  107-117. 

24.  N.  Bleistein  and  R.  Handelsman,  Asymptotic  Expansions  of  Integrals  (Holt,  Rinehart,  and 
Winston,  New  York,  1975),  Ch.  9,  pp.  367-410. 

25.  N.  Bleistein  and  R.  Handelsman,  Asymptotic  Expansions  of  Integrals  (Dover,  New  York,  1986), 
Ch.  9,  pp.  367-410. 

26.  A.  Erdelyi,  “Uniform  Asymptotic  Expansion  of  Integrals,”  in  Analytical  Methods  in  Mathe¬ 
matical  Physics,  R.P.  Gilbert  and  R.G.  Newton,  eds.  (Gordon  and  Breach,  New  York,  1970), 
pp.  149-168. 

27.  R.A.  Handelsman  and  N.  Bleistein,  Arch.  Ration.  Mech.  Anal.  35,  267-283  (1969). 

28.  K.E.  Oughstun,  “Propagation  of  Optical  Pulses  in  Dispersive  Media,”  Doctoral  dissertation, 
The  University  of  Rochester,  Rochester,  N.Y.,  1978,  Ch.  2;  and  K.E.  Oughstun  and  G.C. 
Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics  (Springer- Verlag,  Berlin, 
1994),  Ch.  V,  p.  141. 

29.  K.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994),  Ch.  5,  Secs.  6.5  and  7.2.1. 


Pulsed  Propagation  for  Dispersive  Media 


35 


30.  K.E.  O  ugh  stun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994),  Ch.  5,  p.  168  and  Ch.  7,  p.  281. 

31.  K.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994),  Ch.  6,  pp.  245-248. 

32.  F.W.J.  Olver,  SIAM  Rev.  12(2),  228-247  (1970). 

33.  C.  Chester,  B.  Friedman,  and  F.  Ursell,  Proc.  Cambridge  Philos.  Soc.  53,  599-611  (1957). 

34.  N.  Bleistein,  Commun.  Pure  Appl.  Math.  XIX(4),  353-370  (1966). 

35.  N.  Bleistein,  J.  Math.  Mech.  17(6),  533-559  (1967). 

36.  K.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994),  Ch.  6,  pp.  192-249. 

37.  N.  Bleistein  and  R.  Handelsman,  Asymptotic  Expansions  of  Integrals  (Holt,  Rinehart,  and 
Winston,  New  York,  1975),  Ch.  9,  p.  390;  and  R.A.  Handelsman  and  N.  Bleistein,  Arch. 
Ration.  Mech.  Anal.  35,  267-283  (1969). 

38.  K.E.  Oughstun  and  G.C.  Sherman,  Electromagnetic  Pulse  Propagation  in  Causal  Dielectrics 
(Springer- Verlag,  Berlin,  1994),  Ch.  6,  p.  214. 

39.  L.  Brillouin,  Ann.  Phys.  (Leipzig)  44(10),  203-240  (1914),  Eq.  (16)  on  p.  212. 

40.  V.K.  Rohatgi,  An  Introduction  to  Probability  Theory  and  Mathematical  Statistics  (Wiley  New 
York,  1976),  Ch.  5,  p.  206. 

41.  The  support  of  a  real-valued  function  on  the  real  line  is  the  closure  of  the  set  of  all  points  at 
which  the  function  is  nonzero. 


Appendix 


As  evidenced  in  the  Sections  2  through  4,  the  reflected  and  transmitted  fields  depend  critically 
on  n,  R,  and  T.  This  appendix  presents  the  behavior  of  each  quantity  throughout  the  cut  complex 
plane  (Fig.  2).  Since  the  reflection  and  transmission  coefficients  are  functions  of  the  index  of 
refraction,  n  is  treated  first. 


The  index  of  refraction  is  specified  by  Eqs.  (6),  (7),  (19),  and  (20(a)).  With  u  =  u'  +  iu". 


,  /  m  =  [( +  ‘fz?  +  K  +  *)2]  [(*/  ~  ^)2  +  K  +  *)2]  1  * 
1  1  n  [(«'  +  wj)2  +  (wff  +  «)2][(w'-«j)a  +  K  +  tf)a] 


-7r  +  Arctan(rffi±^ 


—  W3  <  0  and  u"  +  $  <  0 


Arctan 


/  («■>»+<)  \ 
\W-^3)  )  ’ 


7T  + Arctan  (rfeSg 


-f  j  o/  —  u>3  =  0  and  u"  +  S  <  0 

-Arctan  ,  u'  -  u3  >  0  and  u"  +  6  <  0 

0,  g/  —  u>3  >  0  and  u>"  +  6  =  0 

a  =  \  / \  (A2) 

Arctan  ((££&)>  g/  —  «3  >  0  and  w"  +  ^  >  0 

f  >  uf  —  o;3  =  0  and  a/'  +  £  >  0 

a-  -  Arctan  >  <*>'  -  «3  <  0  and  w"  +  £  >  0 

.  it,  u'  -  g>3  <  0  and  u>"  +  6  =  0, 

r  7r  +  Arctan  (^7^)  ,  g/  -  w2  <  0  and  g;"  +  £  <  0 

~2~  j  —  u) 2  ~  0  and  ci/*  H~  ^  <1  0 

2tt  -  Arctan  >  u>  ~  w2  >  0  and  u"  +  £  <  0 

,  _  ^  0,  g/  —  g>2  >  0  and  g/'  +  S  =  0 

Arctan  (ffi-tS)  5  g>'  —  W2  >  0  and  w"  +  £  >  0  1  ' 

y  ?  w*  —  CJ2  =  0  and  uj"  -f*  8  >  0 

7r  -  Arctan  ,  a;'  -  g;2  <  0  and  g>"  +  6  >  0 

u?f  -  v2  <  0  and  a ;"  +  £  =  0. 

One  obtains  the  arguments  a'  and  if)9  by  replacing  u9  —  u?3  and  a;'  —  w 2  with  a?'  +  a>3  and  a?'  +a;2  in 

the  expressions  for  a  and  respectively.  The  representations  for  \n\  and  (a  +  a9  -  ^  -  V>') /2,  given 

by  Eqs.  (A1)-(A3),  are  equivalent  to  the  magnitude  and  argument  of  n  in  Eq.  (6.1.14)  of  Ref.  8. 

In  the  remainder  of  the  Appendix,  the  horizontal  (real)  axis  of  each  plot  is  normalized  with  a 
division  by  u>o;  that  is,  w$  corresponds  to  unity.  Moreover,  all  figures  are  graphed  for  uc  =  1016rad/s 
and  for  Brillouin’s  values  of  6,  and  S .  Consequently,  u>c,  6,  u2^  and  a>3  correspond  to  0.25, 
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1.118,  0.07,  0.9975,  and  1.4984,  respectively.  In  addition,  the  branch  points  (u?H,  w+,  and  u>+) 

are  located  approximately  at  —1.4984  —  0.07i,  —0.9975  —  0.07i,  0.9975  —  0.07i,  and  1.4984  —  0.07 £. 
Although  specific  values  of  the  physical  parameters  are  chosen,  the  analysis  is  independent  of  these 
values,  and  the  figures  typify  the  behaviors  of  n,  R ,  and  T  for  any  choice  of  wCJ  Vo,  6,  and  S  under 
the  constraint  S  <  cjq* 


Even  though  n(w)  is  discontinuous  across  the  branch  cuts,  setting  u;"  =  in  Eq.  (Al)  implies 
that  |n(u;)|  approaches  0  as  u  — ►  and  approaches  +oo  as  w  — ►  Further,  one  can  argue  that 

|n(a;)|  has  well-defined  limits  a s  u>  approaches  each  point  of  the  branch  cuts.  In  particular, 


lim  |n(a;)|  = 

ui— ►a— i6 


(A4) 


for  a  in  [— W3,—  0*2]  or  [<*>2,^3]-  Consequently,  |n(u>)|  can  be  made  continuous  for  u*  in  the  complex 
plane  except  at  u*±,  where  the  limits  exist  but  are  not  finite.  Both  the  3-D  (Fig.  Al(a))  and  cross- 
sectional  (Fig.  A2(a))  depictions  illustrate  the  pronounced  spikes  at  u±  and  the  dips  to  zero  at 
u>±.  In  the  3-D  graph,  the  flat  part  of  each  peak  is  indicative  of  magnitude  clipping  at  the  value  2. 
Figure  Al(a)  provides  a  further  illuminating  glimpse  of  the  behavior  of  |n(u>)|  near  the  branch  cuts 
by  showing  the  cross  section  that  is  slightly  above  the  cuts  at  u"  =  —  0.0699u*o.  Clearly,  away  from 
the  branch  cuts,  |ra(u*)|  is  approximately  unity.  Specifically,  for  |o>|  >  4u*o,  the  calculated  values  of 
|n|  are  within  0.01  of  1. 

In  Figs.  Al(b)  and  A2(b),  plots  of  the  argument  (a  +  a'  -  xf>  -  4>')/2  of  n  indicate  a  jump 
discontinuity  at  each  point  along  the  line  through  the  branch  points  (u"  =  —6).  In  particular,  as 
co"  approaches  -6  from  above,  arg[n(u*)]  approaches  0  for  u'  in  the  intervals  (-00,0*3),  (-0*2, 0*2), 
and  (o>3,+oo),  approaches  -*r/2  for  u>'  in  [-u*3,-a>2],  and  approaches  x/2  for  u'  in  [o*2,^3]-  As  u>" 
approaches  —6  from  below,  arg[n]  approaches  -27r  for  a/  in  (-00,0*3),  (— o*2,o*2),  and  (0*3,  +00), 
approaches  -3ir/2  for  0/  in  [-0*3,  -0*2],  and  approaches  —  5tt/2  for  o*'  in  [o>2,o*3].  These  graphical 
observations  are  consistent  with  the  entries  of  Table  1.  Despite  the  discontinuous  behavior  of  arg[n] 
along  u"  =  —iS,  e’  ar*t71l/2  js  continuous  at  each  point  of  this  line  that  is  not  on  the  cuts.  Of  course, 
one  expects  a  difference  between  the  values  on  the  top  and  bottom  of  the  cuts.  With  the  choice  of 
branch  given  by  Eqs.  (22),  e,argtn]/2  is  -ft  (-t)  for  the  top  of  the  right  (left)  cut  and  is  —i  (-ft)  for 
the  bottom  of  the  right  (left)  cut.  As  u"  moves  away  from  -S,  arg[n]  lies  in  (-f ,  f )  and  approaches 
0  for  u"  >  -6,  and  argfn]  lies  in  (— ^,—  and  approaches  -27T  for  u"  <  —6.  One  consequence 
of  the  preceding  discussion  is  that  cos[(o  +  a'  -  -  i>')/2]  >  0  for  all  u  in  the  complex  plane  and 

is  zero  along  the  branch  cuts.  Geometric  arguments  of  these  observations  are  readily  obtained. 


Since  cos[(a  -f  a'  -  ip  —  $')/2]  >  0,  the  real  part  tir  of  n  is  nonnegative.  In  fact,  n/?  >  0  except 
at  u'±,  and  one  can  rigorously  argue  that  the  limiting  value  is  0  as  u*  approaches  either  branch 
point.  The  plots  for  nR  (Figs.  Al(c)  and  A2(c))  and  |n|  (Figs.  (Al(a)  and  A2(a))  are  very  similar, 
except  that  nR  appears  to  have  steeper  slopes  near  the  infinite  discontinuities  associated  with  the 
other  branch  points  u*±.  As  |u*|  increases  without  bound,  nR  approaches  unity.  For  |u*|  >  4u*o> 
\nR-l\  <  .01.  On  the  other  hand,  tij  is  approximately  zero  for  u  away  from  the  branch  cuts  (Figs. 
Al(d)  and  A2(d)),  and  |n/|  <  .01  for  |u;|  >  4a*o.  Moreover,  the  infinite  discontinuities  of  |n|  at  u>± 
also  appear  in  n/.  Since  Fig.  Al(d)  is  clipped,  the  peaks  in  Fig.  A2(d)  give  an  indication  of  the 
discontinuities. 


Similar  to  Figs.  Al  and  A2,  Figs.  A3  and  A4  exhibit  in  three  dimensions  and  in  two  dimensions 
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four  aspects  (|T|,  arg[T],  Re[T],  Imp1])  of  the  transmission  coefficient.  According  to  Eq.  (65), 

m«)i  =  -7  2  (as} 

y 1  +  |n(«)|2  +  2|n(w)|  cos(arg[n(v)]) 

~lra(a?)l  sin(arg[n(q>)])  1 

1  +  |n(u>)|cos(arg[n(w)])  ’  '  ' 

for  uj  not  along  the  branch  cuts. 

By  Table  1,  arg[n]  takes  on  the  values  =§?,  =f^,  =f,  and  §  along  the  tops  and  bottoms  of 
the  cuts,  which  implies  that  cos(arg[n(w)])  =  0  on  the  cuts.  Further,  |n(o;)|  is  given  by  Eq.  (A4) 
on  the  cuts.  Hence,  |T(w)|  is  well-defined  and  continuous  over  the  complex  plane  with 


arg[T(u>)]  =  arctan 


|2»l  =  2 


along  the  branch  cuts.  Letting  a  — >  ujn  and  a,  — ►  u>3  imply 


lim  l^’(a;)l  =  2  and  lim  |T(u>)|  =  0. 

UJ-+LJ’±  U)— V 


(A7) 


(A8) 


Figures  A3(a)  and  A4(a)  indicate  that  \T\  has  maxima  and  minima  at  the  branch  points  u>'±  and 
oj±,  respectively,  and  quickly  approaches  unity  for  u>  away  from  the  branch  cuts. 


Also  because  cos  (arg[wj)  is  nonegative  in  the  complex  plane  and  is  zero  on  the  cuts,  the 
denominator  in  Eq.  (A6)  is  positive.  Thus  -f  <  arg[T(w)]  <  f ,  and  arg[T(w)]  =  0  when  lies 
along  the  line  oj  =  -id  minus  the  branch  cuts.  Along  the  cuts, 


arg[T(w)] 


Arctan 


j  Arctan 


Top  Right  &  Bottom  Left  Cuts 


Bottom  Right  k  Top  Left  Cuts. 


(A9) 


Clearly,  the  argument  of  T  is  discontinuous  along  each  cut,  except  at  u>  =  u'±  (o  =  ±w3)  where 
arg[T(w)]  =  0.  This  discontinuous  behavior  directly  reflects  the  discontinuities  of  arg[n(o>)]  along 
the  cuts.  Equation  (A9)  also  implies  that  arg[T(w)]  — >  —  j  asw±  along  the  top  right  and 
bottom  left  cuts  and  that  axg[T(w)]  ->  |  along  the  bottom  right  and  top  left  cuts.  The 

7r  jumps  in  the  argument  of  T  at  are  borne  out  by  the  sharp  inner  edges  of  the  two  fin-like 
structures  in  Figs.  A3(b)  and  A4(b). 


From  Eqs.  (65)  and  Figs.  A3(c)  and  A4(c),  one  can  deduce  that  the  real  part  Tr  of  T(u)  is 
continuous  in  the  cut  plane,  that  Tr  can  be  made  continuous  along  the  cuts  by  defining  it  with 
Eq.  (A7)  at  those  points,  that  0  <  Tr(u)  <  2,  and  that  the  maximum  (minimum)  value  2(0)  of  Tr 
occurs  at  u'±  (w±).  Phenomena  that  are  clearly  visible  in  both  figures  include  a  ridge  above  each 
cut,  peaks  about  Wj.,  valleys  about  ,  and  a  leveling  off  to  the  limiting  value  of  1  as  u;  moves 
away  from  the  branch  cuts. 
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On  the  other  hand,  the  presence  of  the  sine  function  makes  T/  discontinuous  along  the  branch 
cuts,  except  at  points  where  the  limiting  value  of  Xj  is  0.  Analytically,  one  can  show  that  |7j(l>)|  <  1 
off  the  cuts  and  that 


Ti(u) 


= 

1  —z&z 


\/(a2  -  W2)(W3  -  “2)> 

\/(a2  “  ul)(ul  ~  fl2)> 


Bottom  Right  h  Top  Left  Cuts 
Top  Right  &  Bottom  Left  Cuts. 


(A10) 


Since  |7/|  <  1  for  to  not  on  the  branch  cuts  and  |T/|  =  1  only  at  a  =  ± [(o>|  -f  u>2 )/2]1/2  along  the  top 
and  bottom  of  each  cut,  clearly  —1  <  T[(u)  <  1  for  all  u.  ± [(<^|  +  —  iS  axe  maxima  along 

the  bottom  right  and  top  left  cuts,  and  ±[(u>3  +  w2)/2]1//2  —  iS  are  minima  along  the  top  right  and 
bottom  left  cuts.  This  behavior  is  apparent  in  Fig.  A4(d),  where  the  dark  curve  approximates  the 
bottom  of  the  cuts  [u/'/cjq  =  -0.0701  or  u"  =  -1.00145]  and  the  light  curve  approximates  the  top 
of  the  cuts  [u>" /u>o  =  —0.0699  or  u"  =  -0.99865].  As  u>  moves  away  from  the  cuts,  T/  approaches 
unity. 


According  to  Eqs.  (20),  n(w)  also  determines  the  behavior  of  the  reflection  coefficient  ( R  = 
Rr  +  iRi).  Figures  (A5)  and  (A6)  provide  3-D  and  cross-sectional  graphs  of  three  of  the  four 
aspects  of  R(u>).  Plots  of  Ri  are  omitted  since  Rj  =  Tj  (see  Figs.  A3(d)  and  A4(d)).  By  Eq.  (A4), 


l*(«)l  = 


1  +  |n(u?)|2  -  2|n(a;)|  cos(arg[n(u>)]) 
1  +  |n(w)|2  +  2|n(w)|  cos(arg[n(w)]) 


(All) 


which  reduces  to  unity  as  uj  approaches  the  top  or  bottom  of  either  cut.  Therefore,  setting  \R\  =  1 
on  the  cuts  insures  that  \R\  is  continuous  over  the  complex  plane.  Since  cos(arg[n(o;)])  =  0,  if  and 
only  if,  u;  is  on  a  side  of  one  of  the  cuts,  \R\  attains  its  maximum  value  of  1  only  on  the  cuts.  As  \u\ 
increases,  |ti|  1  and  arg[n]  — >  0  or  — 2tt  (Fig.  Al(b)  and  Table  1).  Consequently,  the  numerator 
in  Eq.  (All),  hence  |/2|,  approaches  0  through  positive  values.  These  characteristics  are  discernible 
in  Figs.  A5(a)  and  A6(a). 


The  argument  of  R , 


arg[f2(u>)]  =  arctan 


— 2|n(u>)|sin(arg[n(u;)]) 

1  - 1#(«)|* 


(A12) 


has  discontinuities  not  only  along  the  branch  cuts  from  the  discontinuities  of  arg[n(u;)]  but  also  along 
the  imaginary  axis  (u?'  =  0)  stemming  from  the  definition  in  Eq.  (A12).  This  latter  discontinuity  is 
manifested  as  the  wall-like  structure  in  Fig.  A5(b),  which  corresponds  to  a  2i r  jump.  Unfortunately, 
the  jump  discontinuities  along  the  cuts  are  hidden  by  the  surface.  To  perceive  them,  look  at  the 
cross-sectional  curves  on  either  side  of  the  line  u/'  =  —6  (Fig.  A6(b)).  The  light  curve  and  the 
dark  curve  correspond  to  u>,f  =  —0.99865  and  un  =  —1.00145,  respectively.  As  u)  — *  on  the  top 

and  bottom  of  the  right  cut,  arg[fZ(o;)]  — ►  —7 r  and  7r,  respectively;  while  arg[iZ(o;)]  — ►  7r  and  -7 r  as 
u  — *  —  u?4.  on  the  top  and  bottom  of  the  right  cut. 


The  real  part  Rr  of  the  reflection  coefficient  is 


Rr{u) 


_ 1  -  IrcMl2 _ _ 

1  +  |n(w)|2  +  2|n(w)|  cos(arg[n(w)])  ‘ 


(A13) 
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Since  cos(arg[ra(o>)])  >  0,  |.Rr|  <  1,  with  equality  occurring  at  the  branch  points.  Specifically, 
Rr{u>'±)  =  1  and  R(u±)  =  -1.  Moreover,  Rr  has  removable  discontinuities  along  the  cuts;  that  is, 
Rr  is  definable  along  the  cuts  by  taking 

Rr{u)  =  RR(a  -  iS )  =  —  ~  ~  (A14) 

W3  -  U>2 

for  a  in  [-u>3,  -u>2]  or  [u>2,w3].  Thus  the  arcs  in  Figs.  A5(c)  and  A6(c)  connecting  the  values  of  Rr 
at  and  w_  are  segments  of  a  parabola. 


Of  special  interest  is  the  behavior  of  T(w(#))  because  of  its  role  in  determining  the  Sommerfeld 
precursor  As  with  Eq.  (66).  Recall  that  Tr  and  Tj  are  the  factors  in  this  equation  and  are  evaluated 
at  the  distant  saddle  point  ujji,  which  itself  depends  on  0.  Table  A1  provides  values  of  Tr  and 
T/  at  the  u>d\  and  u>m  of  Tables  3  and  4.  The  evaluations  at  the  near  saddle  points  u>ni  are 
included  for  completeness.  As  6  increases  from  1  to  20,  ($))  monotonically  increases,  while 

Ti(u>di(0))  is  a  negligible,  positive  number  for  each  value  of  6.  In  fact,  Tj  never  exceeds  0.014.  In 
the  two  limiting  situations,  $  -+  1+  and  0  -*  +oo,  the  limiting  values  of  T/(wDl(»))  are  1  and  2, 
respectively.  These  limits  are  obtained  from  Eqs.  (65),  Eqs.  (19),  and  Fig.  2  by  noting  that 


cos[(a  +  a'  -  -  t/>')/2]  =  1, 


(A15a) 


sin[(a  +  a'  -  -  V>')/2]  =  0, 


Urn  BE... 

rr 

lim  ~ 
e^+oo  rr 


i, 

o. 


(A15b) 

(A15c) 

(A15d) 
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Table  A1  -  Transmission  Coefficients  Evaluated  at  the  Distant 
(k>r>i(0))  and  Near  (^>^1(0))  Saddle  Points,  Where  0  is  Normal¬ 
ized  Time 


e 

T(um(e)) 

1.000 

1.0000  -  0.00000* 

0.8958  -  0.00000*' 

1.009 

1.0044  -  0.00007*' 

0.8921  -  0.00000* 

1.010 

1.0048  -  0.00008* 

0.8917  -  0.00000* 

1.015 

1.0071  -  0.00014* 

0.8897  -  0.00000* 

1.020 

1.0094  -  0.00020* 

0.8879  -  0.00000* 

1.025 

1.0115  -  0.00027* 

0.8861  -  0.00000* 

1.030 

1.0136  -  0.00034* 

0.8843  -  0.00000* 

1.040 

1.0177  -  0.00048*' 

0.8810  -  0.00000*' 

1.050 

1.0216  -  0.00062* 

0.8779  -  0.00000* 

1.100 

1.0392  -  0.00129* 

0.8643  -  0.00000* 

1.150 

1.0546  -  0.00190* 

0.8531  -  0.00000* 

1.200 

1.0685  -  0.00246* 

0.8434  -  0.00000*' 

1.250 

1.0813  -  0.00296* 

0.8348  -  0.00000* 

1.300 

1.0933  -  0.00343* 

0.8269  -  0.00000*' 

1.400 

1.1154-  0.00429* 

0.8129  -  0.00000* 

1.500 

1.1356  -  0.00505* 

0.8000  -  0.00000* 

1.502735905 

1.1361  -  0.00507* 

0.7994  -  0.00000*' 

1.510 

1.1375  -  0.00512* 

0.7986  -  0.00047* 

1.520 

1.1395  -  0.00519* 

0.7976  -  0.00073* 

1.530 

1.1414  -  0.00526* 

0.7966  -  0.00091* 

1.540 

1.1433  -  0.00533* 

0.7955  -  0.00106* 

1.550 

1.1452  -  0.00540* 

0.7945  -  0.00119* 

1.600 

1.1544  —  0.00574* 

0.7896  -  0.00168* 

1.700 

1.1720  -  0.00638* 

0.7805  -  0.00231* 

1.800 

1.1887  -  0.00696*' 

0.7722  -  0.00276* 

1.900 

1.2046  -  0.00751*' 

0.7646  -0.00311* 

2.000 

1.2197  -  0.00801* 

0.7574  -  0.00340* 

2.500 

1.2872  -  0.01002* 

0.7278  -  0.00435* 

3.000 

1.3441  -0.01140* 

0.7047  -  0.00493* 

4.000 

1.4359  -  0.01284* 

0.6694  -  0.00562* 

5.000 

1.5068  -  0.01323* 

0.6426  -  0.00602* 

10.000 

1.7023  -  0.01094* 

0.5618  -  0.00677* 

20.000 

1.8362  -  0.00688* 

0.4845  -  0.00692*' 

100.000 

1.9647  -  0.00162* 

0.3261  -  0.00591* 

1000.000 

1.9964  -  0.00016* 

0.1689  -  0.00355* 

00 

2.0000  -  0.00000* 

0.0000  -  0.00000*' 
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::: 


Fig.  A1  -  Refractive  index  n  over  the  cut  complex  plane  for  a  singly  resonant  Lorentz  medium  with  o?o  =  4.0  X 
1016s-\  b2  =  20.0  X  1032s-2,  and  S  =  0.28  X  1016s-1:  (a)  |n(w)|,  (b)  arg[ra(u>)],  (c)  Re[n(w)],  and  (d) 
Im[n(u>)].  The  real  (CJf)  and  imaginary  axes  are  normalized  by  Uq. 
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Fig.  A2  -  Cross-sectional  views  of  the  refractive  index  71  for  a  singly  resonant  Lorentz  medium  with  u o  =  4.0  X 

1016s-\  £>2  =  20.0  X  1032s-2,  and  S  =  0.28  x  1016s~1:  (a)  |n(u?)|  along  uff  =  —  0.0699u>o,  (b)  arg[rz(u;)] 
along  u"  =  — 0.0699u>o  (light  curve)  and  a/'  =  — 0.0701a?o  (dark  curve),  (c)  Re[n(u;)]  along  un  =  ~0.0699o;o, 
and  (d)  Im[n(u>)]  along  Un  =  — 0.0699u>o  (light  curve)  and  Un  =  — 0.0701u?o  (dark  curve).  The  real  (a/)  and 

imaginary  (tJu)  axes  are  normalized  by  Uq. 
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Fig.  A3  -  Transmission  coefficient  T  over  the  cut  complex  plane  for  a  singly  resonant  Lorentz  medium  with  u>o  = 
4.0  X  1016s-\  b2  =  20.0  X  1032s-2,  and  6  =  0.28  X  lO16®'1:  (a)  |T(w)|,  (b)  arg[T(«)],  (c)  Re[T(w)],  and 
(d)  Im[T(w)].  The  real  (u>r)  and  imaginary  (u")  axes  are  normalized  by  u-'n. 
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ir«o)l  ar£[r(w)] 


Fig.  A4  -  Cross-sectional  views  of  the  transmission  coefficient  T  for  a  singly  resonant  Lorentz  medium  with  U>q  = 
4.0  X  lO16^1,  b2  =  20.0  X  1032s-2,  and  S  =  0.28  X  lO16^1:  (a)  \T(u)\  along  u"  =  -0.0699^0,  (b) 
arg[T(w)]  along  u)n  —  —0.0699^0  (light  curve)  and  u>n  =  —  0.0701u>o  (dark  curve),  (c)  Re[T(u>)]  along 
Utr  =  — 0.0699u?o,  and  (d)  Im[T(u>)]  along  Un  =  —  0.0699u;o  (light  curve)  and  —  — 0.0701u>o  (dark  curve). 

The  real  (a/)  and  imaginary  (a/7)  axes  are  normalized  by  CJq. 


Fig.  A5  -  Reflection  coefficient  R  over  the  cut  complex  plane  for  a  singly  resonant  Lorentz  medium  with  = 

4.0  X  1016s-\  b2  =  20.0  X  1032s~2,  and  6  =  0.28  X  1016.-1:  (a)  |JS(«)|,  (b)  arg[i2(w)],  (c)  Re[E(o>)],  and 

(d)  Im[i2(o;)].  The  real  (u>')  and  imaginary  (oj")  axes  are  normalized  by  Uq. 
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M*«o)] 


(C) 

Fig.  A6  -  Cross-sectional  views  of  the  reflection  coefficient  R  for  a  singly  resonant  Lorentz  medium  with  Wo  = 
4.0  X  1016s“\  b2  =  20.0  X  1032s~2,  and  <5  =  0.28  X  1016s”*:  (a)  \R(v)\  along  w"  =  -0.0699u;o,  (b) 
arg[i£(u;)]  along  un  —  — 0.0699cJq  (light  curve)  and  Un  =  — 0.0701o;o  (dark  curve),  (c)  Re[jR(w)]  along 
UN  =  — 0.0699u;o,  and  (d)  Im[j?(o;)]  along  Un  =  —  0.0699u;o  (light  curve)  and  UJn  =  — 0.0701tJo  (dark  curve). 

The  real  (a/)  and  imaginary  (u//)  axes  are  normalized  by  U)q. 


